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The issue of non-analytic corrections to the Fermi-liquid behavior is revisited. Previous studies 
have indicated that the corrections to the Fermi-liquid forms of the specific heat and the static 
spin susceptibility scale as T D and T ^ 1 , respectively (with extra logarithms for D = 1,3). In 
addition, the non-uniform spin susceptibility is expected to depend on the bosonic momentum Q in 
a non-analytic way, i.e., as Q D_1 (again with extra logarithms for D = 1, 3). It is shown that these 
non-analytic corrections originate from the universal singularities in the dynamical bosonic response 
functions of a generic Fermi liquid. In contrast to the leading, Fermi-liquid forms which depend on 
the interaction averaged over the Fermi surface, the non-analytic corrections are parameterized by 
only two coupling constants, which are the components of the interaction potential at momentum 
transfers q = and q = 2kF- For 3D systems, a recent result of Belitz, Kirkpatrick and Vojta for 
the spin susceptibility is reproduced and the issue why a non-analytic momentum dependence of the 
non-uniform spin susceptibility (Q 2 lnjQ|) is not paralleled by a non-analyticity in the T— depen- 
dence (T 2 ) is clarified. For the case of a 2D system with a finite-range interaction, explicit forms of 
the corrections to the specific heat (oc T 2 ), uniform (oc T) and non-uniform (oc |Q|) spin suscepti- 
bilities are obtained. It is shown that previous calculations of the temperature dependences of these 
quantities in 2D were incomplete. Some of the results and conclusions of this paper have recently 
been announced in a short communication [A. V. Chubukov and D. L. Maslov, cond-mat/0304381 . 

PACS numbers: 71.10Ay, 71.10 Pm 



I. INTRODUCTION 

The universal features of Fermi liquids and their physical consequences continue to attract the attention of the 
condensed-matter community for almost 50 years after the Fermi-liquid theory was developed byLandau fil. A 
search for stability conditions of a Fermi liquid and deviations from a Fermi- liquid behavior U H 0, B B 0, § , 
particularly near quantum critical points, intensified in recent years mostly due to the non-Fermi-liquid features of 
the normal state of high T c superconductors [9( and other materials. 

In a generic Fermi liquid, the imaginary part of the retarded fermionic self-energy Hn(k,ui) on the mass shell is 
determined solely by fermions in a narrow (~ u>) energy range around the Fermi surface and behaves as p!oLllll[l2lll3| | 

^ = a{^ + {*tY). (i.i) 

Simultaneously, the real part of the self-energy scales as T/ R = But, at small energies (Kramers-Kronig relations relate 
constants A and B via an ultraviolet energy cutoff W ~ Ep). A regular form of the self-energy has a profound effect 
on observable quantities such as the specific heat and static spin and charge susceptibilities, which have the same 
functional dependences as for free fermions, e.g., specific heat C(T) is linear in T, while the susceptibilities Xs(Q,T) 
and Xc(Q,T) both approach constant values at Q = and T = 0. A regular behavior of the fermionic self-energy is 
also in line with a general reasoning that turning on the interaction in D > 1 should not affect drastically the low- 
energy properties of an electronic system, unless special circumstances, e.g., proximity to a quantum phase transition, 
interfere UM 

The subject of this paper is the analysis of the non-analytic, universal corrections to the Fermi-liquid behavior that 
should be present in a generic Fermi liquid. It has been known for some time that the subleading terms in the to— 
and T— expansions of the fermionic self-energy do not form regular, analytic series in uj 2 or T 2 (i.e., w 3 ,oj 5 , etc. for 
£' and w 4 ,o; 6 , etc. for £") 14] . In particular, in D = 3, the power counting shows the first subleading term in the 
(retarded) on-shell self-energy at T — is 16] 

6E R (u) = £ a (u;) - Ej^w) 

= B 3D uj 3 \n(-iuj). (1.2) 
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where B^n is real. For a generic 2 < D < 3, this subleading term behaves as u D . In 2D, it is again logarithmic [TtI 

13 13 s in mm 

i5Efl;(w) = — iB2D<jJ 2 ln(— iw) (1-3) 

where B2D is real. From a formal perspective, the w 2 lnw form of the correction term in 2D implies that at ui — > it 
dominates over a Fermi-liquid, lj 2 -term, i.e., a conventional Fermi- liquid reasoning breaks down. This is true also for 
D < 2, as the correction term scales again as lo d . However, as long as D > 1, KeY, FL (kp, uj) ~ w is asymptotically 
larger at low frequencies than ImE(fcp,w), i.e., fermionic excitations remain well-defined. For a complete set of 
references on this problem see Ref. |24(. 

The singularity in ReS affects directly the subleading term SC(T) in the specific heat C(T) = 7T + SC(T) via ^(j 



*C(T) = 2ttV A 



1 /• d^fc p 0n(w) . 
' " aujuj — Re2j(a;, k)o{u> — e/.) 



T J (2t:)D du 



(1.4) 



where V is the system volume. In 3D, t he p ower counting yields SC(T) oc T 3 lnT [2{|, while in 2D, Re(<5£(w)) oc uj 2 , 
and by power counting SC(T) cx T 2 

Belitz, Kirkpatrick and Vojta (BKV) argued |2^| that the singularity in the fermionic self-energy should also affect 
spin susceptibility an d g ive rise to a singular momentum expansion of the static Xs(Q,T = 0). A similar idea was 
expressed by Misawa (2^ ■ Indeed, the susceptibility is a convolution of the two fermionic Green's functions (a particle- 
hole bubble). For non-interacting fermions, Xs(Q> 0) is given by the Lindhard function which is analytic in Q for small 
Q in all D. The corrections to the Lindhard form are obtained by self-energy and vertex-correction insertions into 
the particle-hole bubble (see Fig. |3J). The diagrams with self-energy insertions can be viewed as convolutions of G 
and Go where G _1 = Gq 1 + S. Substituting the self-energy and expanding in £ and in Q, we obtain 



Sx(Q, 0) = X (Q, 0) - x(0, 0) cx Q 2 I dujde q ^ {UJ ' Q) (1.5) 

(iuj - e q ) b 



Substituting the singular part of S(w) into i|1.5fl and just counting powers, we find Sx(Q,0) oc Q 2 ln|Q| for D = 3, 
and Sx(Q, 0) oc Q ^ 1 for smaller D. (For D = 1, a more accurate estimate yields x(Qi 0) oc In |Q|). 

To verify this reasoning, BKV explicitly computed Sxs(Q,0) in 3D to second order in the interaction, and indeed 
demonstrated |3| that 5x(Q,0) oc <5 2 ln|<5|, in agreement with power counting. Based on this agreement, BKV 
conjectured that power counting should be valid for all D > 1, i.e., the fully renormalized spin susceptibility should 
scale with momenta as Q . 

Another non-analytic behavior was discovered in the analysis of the temperature dependence of the uniform sus- 
ceptibility in 2D. Baranov, Kagan and Marenko (BKM) [30| estimated Xs(Q = 0, T) using a relation between the 
uniform susceptibility and the quasiparticle interaction function [3I3> anc ^ argued that Xs(0,T) is linear in T in 
2D. Chitov and Millis (CM) [3l| later used the same approach, but went beyond estimates and performed a detailed 
analysis of the quasiparticle interaction function and the susceptibility. They also found a linear-in-T dependence. 

Another example of non-analyticity in the leadin g co rrections to a Fermi-liquid behavior is linear-in-T corrections 
to the impurity scattering time in two dimensions |32l |3 0, Hf| . A general treatment of this situation [3(j shows 
that the correction to the residual conductivity of a dirty Fermi liquid depends linearly on the temperature in the 
ballistic regime, i.e., when T is much larger than the level width due to impurity scattering. Unlike the familiar InT 
-dependence of the conductivity in the diffusive regime |37|. this linear T— dependence originates from the singular 
behavior of the response functions of a clean Fermi liquid in 2D. 

Our intension to pursue a further study on singular corrections to the Fermi-liquid behavior is stimulated by 
several factors. First, we want to clarify what actually causes the singularities in the fermionic self-energy, specific 
heat and spin susceptibility. To illustrate the importance of understanding this issue, we note that power counting 
arguments are not rigorous and can lead to incorrect results. Indeed, let's apply power counting to the susceptibility of 
noninteracting fermions, which, we know, is a Lindhard function. Each Green's function of free fermions Go(p-,uj n ) — 
[iuj n — vp(k — kp)]^ 1 scales as one inverse power of momentum and energy (the corresponding dynamical exponent 
zp = 1), so the convolution of the two Green's functions contributes two powers of k — kp in the denominator of the 
integrand for x(Qj0). Expanding up to Q 2 , one then adds two extra powers. The frequency integration eliminates 
one, so there are three powers of momentum left in the denominator. The prefactor for Q 2 should then scale as 

-jr> (1.6) 

where q = p — kp. The lower limit of the integration is of order Q, the upper limit is of order kp. The integral is 
infrared divergent for D < 3 , scales as In \Q\ for D — 3, as IQI 15-3 for 1 < D < 3, and as \Q\~ 2 In |Q| for D = 1. We 
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see that a power counting predicts a singular momentum dependence of the Lindhard function. The true Lindhard 
function obviously does not obey this behavior - it is analytic near Q = for all D. In 3D |38l| . 

XQ {Q,T = Q)= X l D (l-^j (1.7) 

where Xo D = w^p/"" 2 - In 2D, it is just a constant for \Q\ < 2kp [3^.l40| . 

Xo(Q,0)=xt D , Q<2k F , (1.8) 

where Xo° = m/it. In ID 

Xo(Q,0) = xJ D (l + ^), (1.9) 

where Xo D = 2m /irk p. The failure of power counting arguments to reproduce the behavior of the Lindhard function 
clearly calls for understanding under which conditions they do work. The same problem holds also for the self-energy, 
as the singular forms of Eqs. (|1.2fl and (|1.3|) are obtained by power counting, and there is no guarantee that the 
coefficients are nonzero. In fact, CM computed the leading correction to the real part of the self-energy in 2D and 
argued that it vanishes. This would imply that the coefficient Bib in (|1.3(l vanishes, and thus the u 2 lnw in ImE^ is 
absent. Our result will be different (see below) - we will find that B2D is finite. 

Another reason to look more deeply into the physics of singularities is the discrepancy between momentum and 
temperature dependences of the susceptibility. The fact that dynamical exponent zp = 1 would normally imply that a 
non-analytic dependence Sx(Q, T = 0) cx Q D ~ X should be paralleled by a non-analytic dependence of Sx(0, T) oc T D_1 . 
In 3D, this analogy would mean that <5x(0,T) oc T 2 lnT. Misawa |4j| did find a T 2 lnT term in his calculations in 
early 70s. However, later Carneiro and Pethick and recently BKV [2|| argued that the T 2 InT term is actually 

absent in 3D. Several explanations have been put forward to explain this discrepancy. BKV suggested that the absence 
of the T 2 InT dependence in 3D is accidental and should not be regarded as a failure of power counting arguments. 
They conjectured that for a generic D < 3, the T D_1 dependence of x s (0,T) should hold. This conjecture was 
verified numerically by Hirashima and Takahashi [43| for D = 2, but no definite conclusion has been drawn because 
of numerical difficulties. 

As we already said, BKM |3(j and CM [3lJ considered x(0, T) in 2D analytically and argued that the linear-in-T 
term is present. Both groups argued that 8 Xs (0,T) oc T comes from 2kp effects (our results are in full agreement 
with this). BKM also argued that a T— dependence is caused by the singular behavior of the quasiparticle interaction 
function for fermions away from the Fermi surface (in equivalent diagrammatic language - by the singular frequency 
dependence of the particle hole bubble near 2kp). CM found that the linear-in-T behavior is caused not only by 
this effect, but also by the non-analytic temperature behavior of the quasiparticle interaction function for fermions 
at the Fermi surface (in diagrammatic language, by the singular T dependence of the static particle-hole near 2kp). 
The relation between the singularity in the particle hole bubble and non-analyticity of % S (0,T) follows from the fact 
that a generic diagram for for the correction to a Fermi- liquid susceptibility, e.g., diagram 1 in Fig. [21 contains a 
combination 

5 X (0,T)~TJ2 J d 2 kG 3 (k,u; n )Tj2 J d 2 qG (k + q, u n + Q m ) U (q, Q m , T) , (1.10) 

where G(k,oj n ) — (iu) n — e&) 1 is the fermionic propagator. At T — 0, a static particle-hole polarization bubble 

n(g, ui = 0, T = 0) in D = 2 has an asymmetric square root singularity at q — > 2kp + [3^, fiol [ii! Ejj . A 

finite T or finite to soften the singularity and yield H(q,u),T) — H(q, 0, 0) oc -y/max{T, co} in the momentum range 

\q — 2kp\ ~ T/vp [3TI liol |4( | . A simple calculation shows that fermions which contribute to 5xs(0,T) have energies 

of order ~ T and are located in a narrow angular range where the angle between vectors k and q is almost 7r : 
1/2 

7r — 9 ~ (T/Ep) 1 . Using this and assembling the powers, one obtains that <m0, T) oc T. 

In 3D, an analogous reasoning yields the T 2 lnT behavior. CM suggested j3J that previous computations in 3D 
might have missed the crucial 2kp effects and hinted that Misawa may be right in that the T 2 InT term may actually 
be present in 3D. 

In the present communication, we analyze in detail the physical origin of the non-analytic corrections to the Fermi 
liquid and clarify the discrepancy between earlier papers. We obtain explicit results in D = 2 for the fermionic 
self-energy, the effective mass, the specific heat, and for spin and charge susceptibilities at finite Q and T = 0, and at 
finite T and Q — 0. We also verify earlier results for D = 3. 
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We argue that a proper treatment of non-analyticities in the fermionic self-energy and in Xs(Q,0) requires the 
knowledge of the dynamical particle-hole response function. We show explicitly that the non-analyticity in the static 
Lindhard function near 2Uf does not give rise to a non- analytic behavior of the self-energy due to extra cancelations. 
For the spin susceptibility, the computation with the static Lindhard function does yield linear in \ Q\ and T-terms, due 
to 2Uf effects, but with incorrect prefactors. We also demonstrate that non-analytic terms in the self-energy and the 
spin susceptibility can be viewed equivalently as coming either from the non-analyticity in the dynamical particle-hole 
bubble near q = 0, or q — 2fcp, or from the non-analyticity in the dynamical particle-particle bubble near zero total 
momentum. Our results do agree with that of BKV who formally considered only q — contributions. However, we 
show explicitly that they indeed computed all possible non-analytic contributions to the static susceptibility, including 
2kp effects, but just used an unconventional labeling of internal momenta in the diagrams. As an essential step beyond 
the BKV work, we show explicitly that the non-analytic terms in all diagrams for Xs(Q, 0) come exclusively from the 
vertices in which the transferred momentum is either or 2fcp, and simultaneously the total momentum is 0. There 
are only such vertices. They can be viewed as two parts of the scattering amplitude with zero momentum transfer 
and zero total momentum: 

r a ,/3;7,<s(fc, —k; fc, —k) = U(0)S ai Sps ~ U(2kF)8 a s8p 1 . (1-H) 

This restriction to just one scattering amplitude is rather non-trivial, as it implies that non-analytic terms in all 
diagrams for the susceptibility depend only on U(0) and U(2kp) but not on averaged interactions over the Fermi 
surface, as in the BKV analysis. A similar result has been obtained recently for the conductivity in the ballistic 
regime in 2D[36j: for a short-range interaction, the conductivity has a non-analytic T— dependent piece, whose 
prefactor depends only on U(0) and £/(2fcp) rather than on the interaction averaged over the Fermi surface. 

Some of the results and conclusions of this paper have recently been announced in a short communication |47| . 

The paper is organized as follows. In Sec. II we briefly review three known non-analyticities in the response 
functions of a Fermi liquid. In the next four sessions we consider a fermionic system with a contact, i.e., q— independent 
interaction. In Sec. Ill we discuss the leading corrections to the self-energy for interacting fermions in 2D. We show 
that the on-shell self-energy has the form of Eq. (|1 . 3f> with a nonzero B2D , and this gives rise to a linear-in-T correction 
to the effective mass, and to T 2 correction to the specific heat. We show that a correction to the effective mass is 
not observable in a magneto-oscillations experiment due to a peculiar cancelation between two T— dependent terms 
in the self-energy. We also briefly discuss self-energy corrections for D = 3. 

In Sec. IV- VI we consider in detail a non-analytic perturbation theory for the charge- and spin-susceptibilities. We 
use the self-energy calculated in Sec. II along with the dynamical Lindhard functions near q = and q = 2kp and 
the dynamical particle-hole bubble near the zero total momentum as building blocks, and obtain analytic expressions 
for charge- and spin-susceptibilities. More specifically, in Sec. IV we present, for completeness, the expressions for 
the spin susceptibility of noninteracting fermions (the Lindhard function) for various D. In Sec. V we consider the 
susceptibility at T = and finite Q. We present the first analytic calculation of Xs(Q,0) in 2D. We explicitly show 
that it scales as \Q\ and compute the prefactor. These 2D calculations require substantially more effort than in 3D 
since the internal momenta in the diagrams are all of order Q, and one cannot simply expand in Q 2 and then cut the 
infrared divergence of the prefactor by Q, because in 2D the divergence is power law rather than logarithmic. We then 
discuss the 3D case for which we reproduce in a novel way the result of BKV that Sxs(Q, 0) oc Q 2 In |Q|. We explicitly 
verify that non-analytic (|Q|) terms obtained either via a "conventional" approach to treat 2kp contributions and the 
technique invented by BKV are the same. We also discuss briefly the ID case. 

In Sec. VI we consider the static susceptibility at finite T. We show that in 2D, x s (0, T) scales as T with a universal 
prefactor. We also show that the linear-in-T dependence come from two effects: from the thermal smearing of the 
static Lindhard function for particles at the Fermi surface, and from the frequency dependence of the dynamical 
Lindhard function (i.e., from particles outside the Fermi surface). In particular, we show that the linear-in-T piece is 
present in all diagrams for x(0,T), including the ones for which the momentum transfer in the Lindhard function is 
near q — (the linear-in-T terms coming from near q = and q — 2k f are equal). Near q = 0, the static Lindhard 
function is analytic, and a linear-in-T susceptibility comes entirely from the non-analyticity in the dynamical part of 
the Lindhard function. BKM considered only the second source of the 0(T) behavior, CM included both effects. Our 
result differs by a factor of 2 compared to that of CM - we could not detect the reason for the discrepancy. We further 
analyze in detail the physical origin for the linear-in-T term in 2D (and T D ~ X for a general D < 3), and discuss to 
which extent it is related to |Q| D_1 term in x s (Q,0). We show that the physics behind T 1 term in x s (0,T) and 
IQI 13-1 term in Xs(Q,0) is, in fact, different. We discuss how the non-analytic term in x(0,T) T evolves with D, and 
show that for D = 3, Xs(0, T) cx T 2 without an extra logarithmic factor. This agrees with Carneiro and Pethick |42] | 
and BKV results that Xs(0,T) in 3D is free from non-analyticities to order T 2 . We also show that although x s (0,T) 
goes smoothly through D = 2, the 2D case is still somewhat special. Finally, we analyzed charge susceptibility and 
found that non-analytic terms in Xc(Q, T) are all cancelled out, i.e., the first corrections to the Fermi-liquid form for 
the charge susceptibility are all analytic. For a 2D case, this result fully agrees with CM. 
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In Sec. VII we consider the case of a finite-range interaction with q— dependent U(q). We demonstrate that non- 
analytic terms appear in a way similar to anomalies in quantum field theory, and depend only on U(0) and U(2k F ), 
not on the momentum-averaged interaction. We show that at both T = and finite T, the non-analytic correction 
to the self-energy depends on t/ 2 (0) + U 2 (2kp) — U(0)U(2kp), while the total non-analytic correction to Xs depends 
only on U 2 (2kp). We show that the charge susceptibility does not have a non-trivial Q dependence-all non-analytic 
terms from individual diagrams cancel out even when U = U(q). 

In Sec. VIII we present our conclusions. Appendices A-D show details of some calculations. 



II. NON-ANALYTICITIES IN THE BOSONIC RESPONSE FUNCTIONS 



We will demonstrate in this paper that the non-analytic corrections to the Fermi-liquid theory are universally 
related to the Fermi-liquid non-analyticities in the dynamical bosonic response functions. To set the stage, we review 
briefly these non-analyticities. 

There are three physically distinct bosonic non-analyticities in a generic Fermi liquid at T = [H EJ El • The 
first is the non-analyticity in the particle-hole response function, 

n pft (Q,O m ) = - J J d °^\ G(p, u n )G(p + Q, u n + n m ) (2.1) 
at small momentum and frequency transfers. For D = 2, 

i&T°(Q, n m ) = — ( i - '" m| ) ■ (2.2) 

For D = 3, 

ttQ=0, o n mk F /, |Om|, -i vfQ \ 

The zero frequency results: n p ^(0, 0) = m/27r in 2D and n p ^(0,0) = rnkp/^'K 2 in 3D, are the densities of states of 
free fermions per one spin orientation. 

The non-analyticity in the particle-hole bubble at small momenta introduces the dependence of Tl p h(Q — ► 0, u> — > 0) 
on the ratio Vl/vfQ, and eventually gives rise to the emergence of a zero-sound collective mode in a Fermi liquid [l3,[ll|. 

The second is the non-analyticity in the particle-hole response function at momentum transfer near 2fc^. For D = 2 



n^ F (Q,o ro ) = - 







o 


2fc F 







_Q_ 



1/2 x 



where Q = Q — 2k F and 



\ 2vpkp 

<^ 2kp. In the static limit, the non-analyticity is one-sided 



n^(Q,o) 



771 

2^ 



1 



Q - 2k F 



1/2 



for Q > 2k F . 



In D = 3, this non-analyticity is logarithmic and odd with respect to Q 38] . In the static limit 

Q . 4k F \ 



Tj2k F , n ^ TTlkp 



1 - — In ■ - 

, 2fc F |Q| 



(2.4) 



(2.5) 



(2.6) 



The dynamical expression is rather complex in 3D, and we refrain from presenting it. 

The 2k F non-analyticity gives rise to long-range Friedel oscillations of electron density in a Fermi liquid |43 and 
eventually accounts for p— wave pairing in electron systems with short-range repulsive interaction 49] . 
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The third is the logarithmic singularity in the particle-particle response function 

U pp (Q,Q m ) = -[ [ G{ P , u n )G(-p + Q, -w n + Sim) (2.7) 



(2tt 

at small total momentum q and frequency to. In 2D, 



TL PP (Q, fi m ) = — In 5 (2. 



where W ~ £j\ In 3D, the functional form is similar. If the full irreducible interaction between electrons is attractive 
for at least one value of the angular momentum, this singularity gives rise to superconductivity at T = |lflj |. In 
the weak-coupling regime that we will be focusing on, the instability occurs at only exponentially small frequencies, 
and we will neglect it, assuming that the system remains normal down to T = 0. Still, as we will see, a non-analytic 
dependence on the ratio Cl m /i)Fq in n pp (Q, Cl m ) will give rise to a non-analyticity in the self-energy and susceptibility. 

In the rest of the paper we show that these non-analyticities give rise to universal subleading terms in the fermionic 
self-energy, effective mass, specific heat, and static spin susceptibility. 



III. FERMIONIC SELF-ENERGY. EFFECTIVE MASS, SPECIFIC HEAT, AND THE AMPLITUDE OF 

MAGNETO-OSCILLATIONS 

In this Section we obtain non-analytic corrections to the fermionic self-energy and consider how they affect observable 
quantities such as the effective mass and the specific heat. We will mostly focus on D = 2, but for the sake of 
completeness will also discuss the situation in D = 3 and D = 1. We also assume for simplicity that the interaction 
is a contact one, i.e., its Fourier transform is independent of momentum. We will restore the momentum dependence 
of U (q) in Sec. VII. 



A. Self-energy of a generic Fermi liquid 



The (Matsubara) fermionic self-energy is related to the Green's function via 

G~ 1 (k,w n ) = GQ 1 (k,uj n ) + £(fc,w n ), (3.1) 

where Gq (k, u) n ) = %u) n — tk and = {k 2 — klfr) /2m. The two nontrivial second-order diagrams for £(&, w„) are 
presented in Fig. ^ 

For a contact interaction with a coupling constant U, the diagrams a) and b) in Fig. ^yield equal functional forms 
of the self-energy, and only differ in the combinatorial factor resulting from the spin summation and the number of 
closed loops. This factor is equal to (-2) for diagram (a) in Fig.^and to (1) for diagram (b). The result for £(&, oj n ) 
can then be re-expressed as a single diagram Fig^; in which the diamond stands for the interaction vertex iU. In the 
analytic form, we have 

X(k,uj n ) = U 2 G (k 1 )G {k 2 )G {k 3 )8{k 1 + k 2 -k 3 -k). (3.2) 

k 1 ,k 2 ,k 3 

For brevity, we introduced temporarily a "relativistic" notation k = (k, u) n ). The diagram in Fig^j can be equally 
re-expressed either via particle-hole polarization operator H p h(Q, Q TO ), as 

S(fe, w n ) = —TU 2 I j^d G °( k + Q^n + On)IWQ, Cl m ), (3.3) 
or via the particle-particle polarization operator, as 

E(fc, w„) = -TV 2 I Std G o(Q - k, n m - Lu n )Tl pp (Q, Sl m ). (3.4) 

We illustrate the last representation in Fig^i. Here and thereafter uj n — it (2n + 1) T and Q m — l-mnT. 

For definiteness, we will proceed with the form of Eq. I|3.3|) and discuss how the non-analyticity in the particle- 
hole bubble gives rise to the non-analyticity in the fermionic self-energy. To shorten the notations, we will use 
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FIG. 1: a) and b) the two non-trivial second-order diagrams for the self-energy; c) an equivalent form of diagrams a) and b) 
( "sunrise" diagram; d) diagram b) in an explicit particle-particle form. 



H p h(Q, f^m) = n(Q, Sl m ) until otherwise specified. We then show that a non-analytic part of the self-energy can be 
viewed equivalently as coming from the non-analyticity in the particle-particle bubble. 

For the analysis of the specific heat, effective mass and fcrmionic damping, we will need the retarded fermionic 
self-energy Y,ji(k,uj) in real frequencies and at finite temperatures. In some cases, it can be obtained directly from 
Y,{k,uj n ) via a replacement iu) n — ► uj + iS. In general though it is rather difficult to deal with discrete Matsubara 
sums. The approach we adopt here will be to find the imaginary part of the retarded self-energy S^(k, a;). The real 
part of the self-energy £^(fc,a>) is then obtained via the Kramers-Kronig relation. 

Applying the spectral representation 



f(ioj n ) = - f dzM± 

7T J Z — lUJ n 



to l|3.3|) . and using ImG^(k + Q, uj) = —it5(uj — e^q), we find 



E'A(k,w) 



-U 2 



d u Q 
(2nf 



8(n- 



e k+Q )n'A(0,fi) 



, fl , uj + f2 

coth tanh 

IT IT 



(3.5) 



(3.6) 



We first remind a reader how the Fermi-liquid form of S'^(k, u>) is obtained. Suppose that uj <C ep. A simple analysis 
of (|3.6() shows that because of the last term in l|3.6() . typical VI are of order of w, i.e., they are also small compared 
to 6f- The imaginary part of the retarded 11^ (Q, Q) is an odd function of frequency, and hence for small frequencies 
H"(Q, D.) — flF(Q, il). Let's now assume that typical vfQ are much larger than typical f2. Then F(Q, O) w F(Q, 0). 
Substituting this into l|3.6|l . we obtain 



S^(k, W ) 



-IP 



(2n) J 



r 5(e k+Q )F(Q,0) / dnn 



coth ■ 



2T 



tanh ■ 



2T 



(3.7) 



We see that as long as the momentum integral is infrared convergent, it is dominated by large Q ~ hp. The momentum 
integral is then fully separated from the frequency integral and yields a constant prefactor. That typical Q ~ kp 
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also justifies a posteriori the assumption that F(Q, f2) « F(Q,0). The easiest way to do the remaining frequency 
integration is to integrate in a finite range —W < ft < W. Shifting the variable in the second term as fi + tu — ► 0, 
and then setting W = oo we find 

E£(k,w) = C [c 2 + (ttT) 2 ] , (3.8) 

where C is a constant. This is a well-known result in the Fermi-liquid theory |lOj |. 

The form of S^(k, to) given by Ea. (|3.8|) is generic to any Fermi liquid provided that the momentum integral is 
dominated by large momenta Q 3> 0,/vp. Higher order terms in Il'^(Q,il) form a series in £l 2n+1 . If we assume 
that the prefactors depend on Q in a regular way, we obtain higher powers of to 2 and T 2 in E'^ . As we already 
mentioned, this form of Y/^kui) yields, upon Kramers-Kronig transformation, a regular frequency expansion of the 
real part E^(fe, to) = Alu + Blo 3 + where the prefactors are regular functions of T 2 . Of particular importance here 
is the absence of ujT term that would result in a linear-in-T renormalization of the effective mass. It then follows 
that non-analytic corrections to E^ can only emerge if the regular expansion of H'^(Q,fl) breaks down for typical 
momenta that contribute to E'^. This is only possible if H'^(Q,fl) contains non-analytic terms that break a regular 
expansion in odd powers of ft, at least for some momenta. The momentum integration should then show at which 
order of the expansion in ft the prefactor will be divergent enough to make the momentum integral in l|3.7fl infrared 
non-analytic. 

We now show that such non-analytic terms do exist and give rise to non-analytic corrections to the Fermi-liquid 
behavior. One of non-analytic corrections comes from the non-analyticity in Tl(Q, ft) at small Q, another comes from 
the non-analyticity in U(Q,£t) at Q = 2kp. We focus on the 2D case and analyze how these two non-analyticities 
affect the self-energy. We then show that the non-analytic correction to E can be viewed equivalently as coming from 
the non-analyticity in the particle-particle response function. 



B. A non-analytic contribution to the self-energy from Q = 



We begin with the non-analyticity in H"(Q,fl) at small Q. Converting (|2.2|) to real frequencies, we find 

Q 



n£(Q,fi) 



2^ (( VF QY-^yn fOT l°l < v f( 
0, otherwise. 



For notational simplicity we suppress in this subsection the superindex Q — 0. We see that the expansion of 11^, 
holds in powers of Q/vpQ. Obviously, at some order of the expansion, the momentum integral becomes infrared non- 
analytic, which violates the assumption that momentum and frequency integrals in the diagram for the self-energy 
are decoupled. 

In D = 2, this happens already at the leading order in fi. Indeed, substituting (|3.9|l into (|3.6(l . linearizing the 
quasiparticle dispersion as £k+Q = Eft + vpQ cos 6 and integrating first over 9 and then over Q with logarithmic 
accuracy, we obtain 



E£(k,o,) = 



mU 2 
167r 3 i; 



dn ft In 



w 2 



\lu - e k \\2n + (lu - e k )\ 



, CI , uj + fl 

coth tanh 

2T 2T 



(3.9) 



where W ~ Ep is the upper cutoff in the integration over vpq. We see that the momentum integral is infrared-singular 
and introduces an extra logarithmic dependence on frequency. 

The calculation of S^(k, oS) in D = 2 requires certain care as E^(k, w), given by Eq. l|3.9|l . diverges logarithmically 
on the mass shell (ui = However, we will see that this divergence does not affect the real part of the self-energy at 
the mass shell and hence does not affect the specific heat. We therefore proceed in this subsection with the self-energy 
(|3.9|) obtained with the linearized spectrum. In Appendix ^ we consider the mass-shell singularity in more detail 
and show that it is in artifact cured by taking into account either a finite curvature of the electron spectrum or higher 
orders in the expansion in U . 

The frequency integral in (|3.9J) can be evaluated analytically at T = 0. and in the limiting cases at a finite T. At 
T = 0, Eq.lEinj) reduces to (at to > 0) 



E^(k,w) 



mU 2 



dnn In- 



w 2 



\uj - e k \\2Q - (cj - e fe )| 



The integration over f2 is now straightforward and yields 

mU 2 



E£(k, W ) 



l6ir 3 v F 



l ( Ml W 

— (U — 6k) } m -. r 

4 V J \uj + e k \ 



- {lu 2 + ^{uu-e k f y\:i 



W 



(3.10) 



(3.11) 
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where the . . . represent the regular lu 2 term. Away from a near vicinity of ui 
irrelevant (to logarithmic accuracy) and T,"(k,uj) can be written as 

E£(k,w) = £?(k,w) + £2(k,w), 
mU 2 , , W 



E?(k, W ) = 



167T 3 Z)|, 

mil 2 



J 1 In 



|w + e fe | : 



-6^, the term with [us — e^) 2 is 

(3.12a) 
(3.12b) 



167T 3 w|, 



2 In- 



w 



\u-e k \' 

We see from (|3.12al that for e k ~ w, both terms scale as cj 2 lnw. In particular, at e/. = 0, 



Li) k,u = In ■ 



(3.12c) 



(3.13) 



Tracing Eq. (|3.12a(l back to (|3.9|1 . we observe that the first term S"(k, u>) comes from the il-dependent part of the 
logarithm in (|3.9|l . and the second term S' 2 '(k, lu) comes from the Jl-independent part of the logarithm. We see that 
for £ 2 '(k, to), the factorization of the momentum and frequency integrations still holds, and as in a Fermi liquid, 
the momentum integral just adds an overall factor that logarithmically depends on the external lu and e k . On the 
contrary, for E"(k,w), the momentum and frequency integrals are coupled. 

The zero-temperature result for the self-energy can be also obtained directly in Matsubara frequencies, without doing 
the analytic continuation first. Expanding in small momentum transfer Q, we have for the Matsubara self-energy at 
T = ' 



E(fc,cj n )| T =o 



mil 2 
' 8tt 4 



W/v F 



QdQ 



(10 



1 



o vfQcosO + e k - i(u n + Q m ) 



i a, 



y/(v F Qy + ni 

The integration over 6 is elementary and yields 

TflU 

S(fc,w„)| T=0 = 



'8n 3 v 2 F 



|sgn(w„ + 



dx- 



y/x 2 + fl 2 n y/x 2 + (tU n + Q m + ie k f 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



where we introduce x — vpQ- Performing finally the integration over x, we obtain with logarithmic accuracy, for 
w n > 0, 



E(fc, (x>„)|t=o 



mU 



16ir 3 Vp 



. mU 2 
bj 2 n + - (u) n + ie k ) z ) In 



W 



u> n + ie k 
W 



In- 



W 



2Q, r 



uj„. + ie k 



+ w4- 7 (w + ie fc ) > 



W 



(3.18) 



Continuing to real frequencies, (ioj n — > w + iO), we indeed obtain H3.11|l for S^.. The Matsubara self-energy can also 
be partitioned into T,i(k,Lu n ) and £2^, u> n ). The first term is singular near the mass surface, while for the second we 
have (to logarithmic accuracy) for a generic e k /u) n , 



v (u \\ .mU 2 
S 2 (k, uj n )\ T=0 = „ „ to. 



W 



2 In — 



Continuing to real frequencies, we obtain 

E 2 (fc, u))\ T =o 

At finite T, instead of (|3~TU|) we have 

mU 2 



16tt 3 4 ™ u. 



mU 2 a , 7T , W\ 

— -3-2-w --sgnw + zln— . 



E^(k,w) 



167T 3 W|, 



O In 



W 2 



\u> - e fc ||20 + (lu - e k )\ 



n uj + Q 

coth — - — tanh — — — 
2T 2T 



(3.19) 



(3.20) 



(3.21) 
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a) b) 

FIG. 2: a) q = contribution to the self-energy; b) q — 2kp contribution to the self-energy. 

It is again convenient to split the self-energy into two parts, £"((*;) and X' 2 '(w), coming from Independent and Cl- 
independent pieces of the logarithm in (|3.21|) . For the O- independent part of the logarithm, the frequency integration 
is the same as in a Fermi liquid, hence 



S'/(k,o;) 



[w + (7rT) J In 



to - ek\ 



For the second term, we have 



mU 
16n 3 v F 



dQfl In - 



ir 



|2n+(w-e fc )l 



, tt , w + fi 

coth tanh 

IT 2T 



(3.22) 



(3.23) 



In this last term, the dependence on the ratio u)/e% is not singular and can be neglected, to logarithmic accuracy. 
Using series representations for the hyperbolic functions we can then re-express the r.h.s. of l|3.23|l as 



4V) = -^([<^W]in ( r/i)W,(^)) 



(3.24) 



where A is a constant, and 



/(*) = 0.79 + V [ dy tanh ^ (yln -g— + - - In ±±\ 
J 2 V \V !| V 12/ — 1 1 



(3.25) 



One can easily make sure that the expansion of f(x) holds in even powers of x. At large x, f(x) w lnx, i.e., at cj ^> T, 
this expression reproduces oc a; 2 ^^. At small x, i.e., at cj <C T, f(x) 0.79 + 0.35a; 2 . 



C. A non-analytic contribution to the self-energy from q m 2kF 



We next consider a singular contribution to S^(k, u>) from momentum transfers close to 2kp. To perform compu- 
tations along the same lines as for Q near 0, we would need to know the form of II(Q, f2) at finite and T, which is 
rather involved. However, we actually would not need this form at all, as we demonstrate that the contribution to 
the self-energy from Q m 2kp is exactly the same as S2(k, ui) defined in (|3.12a|l . The most straightforward way to see 
this is to go back to a diagram representation of the self-energy in terms of three fermionic propagators (Figlfl;). In 
analytical form, the "q — 0" contribution to the self-energy is 



E 9=0 (/c) = U 2 



jD + l, 



(2tt 



D + l 



d D+1 



(27T) 



-D+l 



Gk+qGpGp+q- 



(3.26) 



where q is assumed to be small. We again use "relativistic" notation k = (k,w) and q = (Q,f2). Integrating over p 
first, we obtain 



s9 =0 = _ V 2 



d D+1 c 



-G k+q Il(q), 



(3.27) 



where 11(g) is a particle-hole bubble at small momentum and frequency. This expression we used in the previous 
subsection. We found that two singular contributions to S 9=0 , Si(k, u) and E 2 (k, w), and that S2(k, u>) comes from 
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the momentum region where two of the internal momenta are close to — k and the third one is close to k, i.e., from 
the range of p which are nearly antiparallel to k. Since p + k are small (of order of external momenta), we can relabel 
the momenta as shown in Fig. |2b and re-express £2(k, ui) as 

f d D+1 q f d D+1 q' 

T, 2 (k) = U 2 / D+1 / - D+1 G k+q G- k+ql G- k+q+q , . (3.28) 

J (27T) J (27T) 

where now both q and q' are small. Integrating over q' first, we obtain a conventional expression for £2^) i n terms of 
the polarization operator with small momentum transfer. On the other hand, changing the order of integration and 
integrating over q first, we obtain 

E 2 (fc) = -U 2 J G- k+q ,fl(2k - q>) (3.29) 

where 

fl(2fc -q') = - [ d ° + ' q +1 G k+q G- k+q+q , . (3.30) 

J (27T) 

In general, II(2fc — q) is not equivalent to the polarization bubble Tl(q) with momentum near 2kp, as our re-writing is 
only valid if internal q are small. However, the singular parts of the two bubbles coincide because the singular part in 
n(Q ~ 2kp, SI) (proportional to \J\Q — 2kp\9(\Q — 2kp\) in the static case) comes from the momentum range where 
the two internal momenta in the particle-hole bubble are close to ±k, i.e., from exactly the same range that is covered 
in II(2fc — q'). We show this explicitly in the Appendix iBl This equivalence implies that the r.h.s. of (|3.29|) is just 
the singular part of the "2/cf" contribution to the self-energy. We see therefore that £ 9=2fcF (fc) = Y,2~°(k). The total 
self-energy is then 

£(fc) = £ 9=0 (fc) + Y, q=2kF = Si(fc) + 2£ 2 (fc). (3.31) 

For momentum-dependent interaction U = U(q), the computation of the 2kp— contribution requires more care and 
we present it in Sec. VII. 

That the 2fcj?-singularity comes from nearly antiparallel internal fermionic momenta has been implicitly used in 
the Kohn-Luttinger analysis of superconducting instability with large angular momenta of Cooper pairs [48j. In the 
context of corrections to the Fermi-liquid theory, Belitz et al. [2^] argued that all singular contributions to the spin 
susceptibility can be described as small q effects, although they did not emphasize that some of their small q effects 
are in fact equivalent to 2kp contributions in conventional notations. 

That both q = and q = 2kp singularities in the polarization bubble contribute to the self-energy was first 
emphasized by CM [3l]]. However, the relative sign of the two terms is different in their and our calculations. We 
found that the singular terms add, while they argued that singular contributions from q = and q = 2kp cancel each 
other. Since the interplay between q — and q — 2kp contributions to the self-energy is crucial to the issue of whether 
or not there is a T 2 -term in the specific heat and linear-in-T term in the effective mass (CM argued that both are 
absent due to cancelation between q = and q = 2k p terms), we present in Appendix B|Clan explicit computation 
of the 2kp contribution to the second-order self-energy at T = 0. This calculation confirms that Yfl =2kF = E 9=0 . 



D. An alternative analysis, in terms of H PP (Q,Q,) 



We next demonstrate that the backscattering non-analyticity in the fermionic self-energy can be viewed equivalently 
as coming from the non-analyticity in the particle-particle bubble at small total momentum and frequency. This readily 
follows from our consideration of the "2/cf" diagram. Indeed, since both q and q' are small, the full self-energy can 
be re-expressed as 

r d D+l n r d D+1 a' 

Performing the same analysis as in the previous section, we observe that the deviation from the Fermi-liquid form 
of E is only possible if the expansion of Hp p (Q,£l) in odd powers of SI breaks down due to infrared divergences 
of momentum dependent prefactors. This is precisely what happens in H pp (Q,tt) given by l|2.8|l as the frequency 
expansion holds in fl/vpQ, i.e., the prefactors are non-analytic at vanishing Q. We emphasize that the logarithmic 
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divergence of n pp (Q, 0) at vanishing Q and is by itself not essential; what matters is a non-analytic dependence on 
the ratio Q/vpQ. 

We see, therefore, that the non-analytic piece in the self-energy can be viewed equivalently as coming from a non- 
analyticity in the particle-hole bubble, or from a non-analyticity in the particle-particle bubble. To further verify this, 
we explicitly compute in Appendix lOl the non-analytic part of £(fc) at T — using the "particle-particle formalism" 
and indeed find it to be equal to the non-analytic 2E2(fc,w) that we obtained in the "particle- hole formalism," i.e., 

£W=°)(A0 = 2£ 2 (fc). (3.33) 

The term Si(fc,w) can be also reproduced in the particle-particle formalism, but this contribution comes from large 
q + q' w 2k, and we refrain from re-deriving this piece. 

Our results on this issue again disagree with those by CM |2jJ. They performed a complimentary analysis of the 
self-energy based on the evaluation of an effective vertex function to second order in U , and argued that there is 
a cancelation between non-analytic contributions coming from the 2kp non-analyticity in the particle-hole channel 
and the non-analyticity in the particle-particle channel. We, on the contrary, find that the contribution from the 
particle-particle non-analyticity is twice the "2/cf" contribution from the particle hole channel. 

Summarizing the results of the last two subsections, we see that the non-analytic part of the fermionic self-energy 
in 2D consists of two parts. The first part, E"(fc), comes from forward scattering. It has the same functional form, 
uj 2 + (irT) 2 , as in a Fermi liquid, but the prefactor logarithmically depends on u> — e&. The second part, E 2 '(fc), comes 
from the processes which involve the scattering amplitude with near-zero total and transferred momentum. This 
E' 2 '(fc) has a non- Fermi- liquid form, and can be equally attributed to the Q — non-analyticity in the particle- hole 
polarization bubble, or to the 2kp non-analyticity in the same bubble, or to the Q — singularity in the particle- 
particle bubble. In the next section we show that only £2^) actually contributes to the thermodynamics. 



E. Effective Mass and Specific Heat 



We first use the result for E" obtained in Sec lIII Bl and compute the real part of the self-energy on the mass shell. 
We then use E'(tj = e k ) to find the effective mass and specific heat. 
The Kramers-Kronig relation on the mass shell is 



dE 



Y 1 "{E,e k =uj) 
E-uj 



We begin with £i(fc). Substituting E"(k, uj) from (|3.12b|) into (|3.34(l . we find that on the mass shell 



Ei(fc,cj)| w=efc = 



16tt 4 v 2 f 



V 



dz 



z 2 + (ttT) 2 



ln- 



TF 



(3.34) 



(3.35) 



By dimensional analysis, the integral in l|3.35|l is of order us 2 . However, the prefactor in front of u> 2 turns out to be 
zero. The easiest way to see this is to evaluate the integral in finite limits — W < z < W and to search for the universal 
term that would be independent of W. Performing elementary manipulations, we find that £i(w) does not contain 
such a term. Foreshadowing, we note that the same result holds for the static spin susceptibility which we discuss in 
detail in Sections II V Al and II V Bl We will see there that the inclusion of the T<2(k,uj) into a particle- hole bubble with 
external momentum Q yields a non-analytic \Q\ term in Xs{Q)- On the contrary, the susceptibility diagram with an 
extra Ei(k, uj) scales, in Matsubara frequencies, as 



Sxoc duj n uj n / de k 



In [W/ (e fc -tw n )] 



(e fe - iuj n ) 2 {(e k - iuj n ) 2 - (v F Q) 2 



(3.36) 



By power counting, the leading Q dependence of the integral should be \Q\. However, a straightforward computation 
shows that the prefactor again vanishes. The outcome of this analysis is that the divergence of E"(k, a;) on the mass 
shell does not give rise to non- analytic corrections to Fcrmi-liquid form of the thermodynamic observables. 

We next consider E 2 (k, u>). Substituting E 2 '(k, u>) from Ea. H^.21|) into Ea. lj3.34jl . we obtain after simple manipula- 
tions 



E 2 H = 



mU 2 
16tt 4 v 2 f 



dnn v 



dE 



E 2 



, n , n+E 

coth tanh 

2T 2T 



LU 



2tt + E-uj 



2n + E 



- uJ 



hi 



\{2Q + E) 2 -uj 2 
W 2 



(3.37) 
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Integrations over f2 and E can be performed exactly. We give the details of this calculation in Appendix[D]and present 
just the results here. At T — 0, we obtain 



E 2 H 



mU 
'32ir 2 v 2 F 



u)\uj\. 



(3.38) 



This coincides with Eq. (|3.20[1 obtained via analytic continuation of the Matsubara self-energy. 
In the opposite limit of small uj/T, we have 



s' 2 H= - 



8tt 2 ? 



ujT. 



(3.39) 



As the self-energy in this region is linear in u>, Eq. 13.39|l implies that the effective mass of subthermal quasiparticles, 
i.e., with ui «T, scales linearly with T. Using the fact that the full £(fc,w) = £i(fc,CL>) + 2S2(fc,w) and that £i(fc,w) 
does not contribute to thermodynamics, we obtain 



to* (T) = to* (T = 0) 1 -21n2 



mU 



(3.40) 



This result disagrees with CM-they argued that the linear- in-T term in the mass renormalization is absent. 

In a very recent study Das Sarma, Galitski, and Zhang |50j did find a linear-in-T correction to the effective mass 
for the Coulomb interaction in D = 2. Although the sign of their linear-in-T term is opposite to that in Ea. (|3.40J) . 
we believe that there is no contradiction here as there are no general restrictions on the sign of the prefactor. It is 
therefore quite possible that the sign of the 0(T) term is different for short- and long-range interactions. Note in 
this regard that the effect of the interaction on the effective mass is different for these two cases even at T — : a 
short-range interaction increases to*, while the Coulomb interaction decreases to* in the limit r s «l [lflj . 

For generic u>/T, the non-analytic part of the full S'(w) can be cast into the following scaling form 



where 



mU 
10 2 2 U M9 I 



r^2 



12 



Li, 



(3.41) 



(3.42) 



and Li2(x) is a polylogarithmic function. 

Note that g{oo) — 1 and g(x « 1) Ri 41n2/x. Substituting these limiting expressions into Ea. H3.41[) we indeed 
reproduce Ea. (|3~3%j) and Eq. (|33^l . 

The full functional form of g(x) is required for the computation of the specific heat, as the frequency integral for 
C(T) given by Ea. (|1.4fl is confined to lu ~ T. Previous work [2(| on C(T) used only the T = form of the self-energy 
and hence yielded incorrect prefactors. Substituting our result for S' into l|1.4l) we obtain in 2D 



SC(T) = C FL 



A8K ( mU 



4n 



T 



where Cfl{T) = mnT/3 is the Fermi gas result for the specific heat and 



K = 



dxx 
cosh 2 x 



h Li 2 f-e 

12 V 



= 1.803. 



(3.43) 



(3.44) 



As it was to be anticipated, the non-analytic correction to the fermionic self-energy gives rise to the T 2 -term in 
the specific heat. It is essential that this non-analytic term comes only from fermions in a near vicinity of the Fermi 
surface and is thus model-independent. The same is true for the linear-in-T correction to the effective mass. In other 
words, the leading corrections to the Fermi- liquid forms of to and C(T) are fully universal. 

The T 2 -dependence of the correction to the specific heat agrees with the results by Coffey and Bedell [27| and 
Misawa j2(|. However, Coffey and Bedell did not explicitly compute the prefactor and apparently only included small 
momentum transfers (i.e., no 2kp effects). Misawa did compute the prefactor, but he neglected the temperature 
dependence of the fermionic self-energy. We found above that this T dependence cannot be neglected, and our 
prefactor disagrees with that by Misawa. 
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F. Amplitude of quantum magneto-oscillations 



In previous sections, we found the general form of non-analytic corrections to the real and imaginary parts of 
the self-energy. We now discuss whether these corrections can be observed experimentally via magneto-oscillations. 
Naively speaking, one might have expected the finite quasiparticle relaxation rate, T 2 lnT, to damp the amplitude 
of the oscillations as a contribution to the "Dingle temperature", whereas the T— dependent effective mass might 
affect the thermal smearing factor. However, we argue below that quadratic and quadratic-times-log terms in the 
self-energy are not detectable by measuring the amplitude of magneto-oscillations in D = 2. 

In the Luttinger formalism [5J|, the amplitude of the k th — harmonic of magneto-oscillations is given by 



A k = — ^ — ex P 

w„>0 



2nk [uj n — i£ (uj n , T)] 



(3.45) 



where f2 c is the cyclotron frequency. It is essential for our consideration that the amplitude is determined by the self- 
energy in the Matsubara representation rather than by the real and imaginary parts of the retarded self-energy [55| . 
By itself, £^, and £^, determine the fermion dispersion and lifetime, respectively; however in l|3.45|l this distinction is 
lost. 

The assumption made in deriving (|3.45(1 is that the dependence of the self-energy on the magnetic field can be 
neglected. In 3D, this assumption is well justified as the effect of the magnetic field on the self-energy yields corrections 
to Ak which are small in , where N — €f/Q c ^ 1 is the total number of Landau levels. In 2D, however, the effect 

of the magnetic field is non-perturbative, and at T — and in the absence of disorder, the field-induced oscillations 
of the self-energy are as important as the oscillations of the thermodynamic potential itself [56|. Ea. l|3.45(l is then 
only applicable as long as oscillations of the thermodynamic potential are exponentially small due to either finite 
temperature and/or disorder. In this paper we disregard effects of disorder (considered recently in [5?|l. thus the 
amplitude is only controlled by the finite temperature. In this case, the restriction of the small amplitude in its turn 
implies that the sum over Matsubara frequencies in (|3.45|) can be truncated to only the n = term. Notice that 
this restriction is mandatory in D = 2 within the Luttinger formalism but depends on the choice of experimental 
conditions in D = 3. The amplitude of the first (largest) harmonic then simplifies to 



4tt 2 T 

Al = ^7 exp 



2tt [ttT - ££ (ttT, T)] 



(3.46) 



The temperature enters the Matsubara self-energy £(w n ,T) in two ways: first, as the Matsubara frequency, and 
second, as the physical temperature determining the thermal distribution of the degrees of freedom. For the lowest 
frequency, ujq — TtT, the interplay between the two effects leads to a peculiar cancelation. 
Indeed, consider for a moment a generic Fermi liquid, for which 



S(w n , T) 



1 ) iuj n + iC 



(3.47) 



where C is a constant, . . . stand for the higher order terms [O (e 3 , T" 3 )] , and m* /m has a regular expansion in powers 
of T 2 . The analytic continuation of (|3.47() to real frequencies yields the correct retarded self-energy (Q) . We see that 
the second term £ (w„,T) vanishes for u> n = ±7rT, i.e., the self-energy that enters into the formula for contains 
terms only of order T 3 and higher. In other words, the quadratic in T piece present in the imaginary part of the 
retarded self-energy and associated observables, does not affect the amplitude of magneto-oscillations, which to order 
T 3 is given by 



Ax 



4n 2 T 



■ exp 



2tt 2 T 
fit 



o* = HLfi c 

m* 



(3.48) 



where m*/m is a regular mass renormalization which comes from fermions far away from the Fermi surface. This 
rather remarkable result was previously obtained specifically for electron-phonon interaction and is known as a "Fowler- 
Prange theorem" jE^ . 

We found that a similar cancelation occurs also for our self-energy in D = 2. To logarithmic accuracy, the second 
term in Ea. H3.47JI is replaced by 



£ (w„,T) 



-iCT ^sgn( 



flm.) \&rn.\ In ■ 



i a. 



w 



(3.49) 
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where C is a real constant, and the factor of sgn(w„ — fi TO ) resulted from the angular integration of the Green's 
function. A simple transformation of the Matsubara sum reduces £ (cj n , T) to 

n,„ =o 

Expression 1|3.5L)[1 obviously vanishes for uj n = irT, i.e., therefore S (ttT,T) in 1)3. 46[) does not contain a contribution 
from £. Due to this cancelation, the exponential factor in A\ does not contain terms of order T 2 lnT. A more 
detailed analysis |57j| . shows that T 2 terms are also absent, i.e., both quadratic terms and quadratic-times-log terms 
in the self-energy (and thus the linear- in-T effective mass [Ea. (|3.40ll ]) are not observable in a magneto-oscillation 
experiments. 



IV. SPIN AND CHARGE SUSCEPTIBILITIES 



We next proceed to the analysis of the corrections to the Fermi-liquid forms of spin and charge susceptibilities. 
The charge and spin operators arc bilinear combinations of fcrmions: 



C(q) = 53 4+q, Q c k,. 



(4.1) 



for charge, and 



S(q) = ^ ff, 



/3 c k+q,a C k,/3 



(4.2) 



k,a,/3 



for spin. The corresponding susceptibilities for a system of interacting fermions are given by fully renormalized 
particle-hole bubbles with side vertices 



a,f3i 



(4.3) 



where c and s refer to charge and spin, respectively. 

For non-interacting fermions, the spin and charge susceptibilities are equal and given by the Lindhard function that 
coincides, up to an overall factor, with the polarization operator Tl(Q, O m ): 



x c o(Q,n m ) = x5(Q,fim) = 2n(Q,fi m ), 

where Xo(Q, ^m) = [Xo(Q, ^m)]u and i = 1,2, 3, and 

d D k 



(4.4) 



n(Q,n m ) = -T53 

m 

G (k,w„) G (k + Q 



(4.5) 



At T = 0,the charge and spin susceptibilities can be evaluated exactly for any Q and fl m . In the static limit, 
= 0, they acquire particularly simple forms. For D = 3, we have |38| 



3D 



1 4fc F - Q 2 Q + 2k F 

2 + 8Qk F U \Q-2k F \ 



X c o (Q,0) = Xo(Q,0) = xt 

where Xo D = nikp/n 2 . In D = 2, the corresponding expression is |39ll4C 

Xo(O.O) = xg(Q,0)=x^, Q<2fc F ; 



(4.6) 



21) 



1-1 



44 x 1/2 * 

n2 



where Xo D — m l^- In ID, we have [5E 



xS(Q,0) = xg(Q,0)= X ^ In 



k F 



k F 



(4.7) 



(4.8) 
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where Xo D = 2/(7TUf). As it was mentioned in the Introduction, Xg S (Q, 0) is analytic in Q for small Q in all dimensions. 

The first nontrivial corrections to Xo' S (Q;0) come from the diagrams presented in FigOU These diagrams represent 
self-energy and vertex-correction insertions into the bare particle- hole bubble |2^. Diagrams 1 — 5 are nonzero for 
both Xs and Xc- Diagrams 6 and 7 are finite for Xc, but vanish for Xs upon the spin summation (X) Q a aa = 0)- 
The internal parts of all diagrams contain fermionic bubbles: particle-hole bubbles for diagrams 1, 2, 3, 5 and particle- 
particle bubble for the diagram 4. 

In the next two sections we analyze the form of the static susceptibility first at a finite Q and zero temperature, 
and then at finite T and Q = 0. 

A. Spin and charge susceptibilities at finite Q and T = 

As in Sec. Ill, we assume that the interaction is independent of momentum. We explicitly computed all 7 diagrams 
Fig. Inland found that each of the diagrams (except for diagrams 6 and 7 which vanish identically for the spin channel) 
contributes a correction Sx(Q,0) oc \Q\, and that this non-analyticity is a direct consequence of the dynamical 
singularities in the particle-hole and particle-particle bubbles. 

1. D = 2 

As we mentioned in the Introduction, the calculation in D = 2 is more difficult to perform than in D — 3 because 
all typical internal momenta and energies are of the same order as the external ones (Q and vpQ, respectively); thus 
no expansion is possible. In 3D, where Sx(Q, 0) oc Q 2 InQ, typical internal momenta are larger than external Q, and 
one could expand the integrand in Q 2 and evaluate the prefactor to logarithmic accuracy. 

We begin with diagram 1 which represents the self-energy insertion into the particle-hole bubble. This diagram 
yields the same contribution for spin and charge channels, so we will drop the subscript and denote xi = Xis — Xic- 

An analytic form of diagram 1 in the Matsubara representation is given by 

S X i(Q,Q) = -8U 2 J Gg(k, cj)G (k + Q, u,)G (k + q, u + fl)II( g , O). (4.9) 

The combinatorial factor of 8 includes two factors of 2 due to spin summation and an extra factor of 2 associated with 
the fact that the self-energy can be added to any of the two fermionic lines in the bubble. Non-analytic contributions 
to Sxi(Q, 0) come from two regions of momentum transfers: q near zero and q near 2kp. Since we have already shown 
in Sec. IIIII that the contributions to the self-energy from these two regions are equal for a contact interaction (up to a 
forward scattering piece in S 9=0 that, as we demonstrated, does not contribute to \Q\ term in the susceptibility), we 
do not have to calculate the q — and q = 2/cp contributions to %i(Q, 0) separately-the two are just equal: 

^r (Q,o) = ^r 2fcF (Q,o). (4.io) 

This implies that we only have to compute Sx1~°{Q,0), the full Sxi(Q, 0) will be twice that value. To be on a safe 
side, we verified this reasoning by explicitly computing 5x\~ 2kF {Q, 0). We present the calculations in Appendix lEl 
We indeed found it to be equal to Sxl~ (Q, 0). 

We now compute SxliQ, 0) Since the non-analyticity in xi(Q> 0) is expected to come from the vicinity of the Fermi 
surface, the fermionic spectra e^, tk+q and fifc+Q can be expanded to first order in k — kp\ 

£k = vp(k - hp), £k+Q = £fe + vfQ cosOi, e^+q = e k + t'Fq cos 9 2 - (4-H) 
Substituting this expansion into Ea. l|4.9|l and performing elementary integrations over k, oj, and 6\, we obtain 

«xr°(o»o) 



2mU 2 f°° f°° 
— / qdq / ndnU(q,fl) (4.12) 

I" JO JO 

o (iil -v F q cos 6 2 ) 2 J(v F Q) 2 + (ft + iv F q cos 6 2 ) 2 ' 



where Tl(q, lo) at small q and u> is given by Ea. H2.2|l . Rescaling the remaining variables as q = q/Q, & = £l/(vpQ) 




FIG. 3: Each of the seven diagrams in this figure give singular corrections to spin and charge susceptibilities. 
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and introducing polar coordinates as q — r cos 0, ft = r sin 0, we obtain from i|4.13|l 

?mTl 2 \<l\ f*/ 2 

5x\ (Q,0) = / d0sin0cos0n(0) (4.14) 

"F JO 

/" 1 1 

d02 / 



(cos cos 6» 2 - i sin 0) 2 ^/l + r 2 (sin + i cos cos 6» 2 ) 2 ' 

where 11(0) = (m/27r)(l — sin0). The upper limit of the integral over r is r max = 0(kp/Q) S> 1. The integration 
over r is straightforward and yields 

*Xr°(Q.O) = ^T^ sin cos 011(0) r^ 2 - i — — (4.15) 

7T 4 Uf Jo Jo (COS0 COS P2 — 1 sin 0J 

1 max 

) 2 (sin0 + icoscf) cosfo) 2 — \Q\ 

As Q r m ax ~ the dominant piece in ^Xi _0 (Q,0) comes from high energies and accounts for the non-universal 
correction to the uniform susceptibility x(0, 0). We, however, are interested in the first subleading term which scales 
as \Q\ and does not depend on r max . Performing the integration over 02, we obtain for this universal contribution 

<5xH°(<3,0) = -H^M r /2 d0 sm 2 0cos0 (5 sin 2 0-3) 11(0). (4.16) 

7T Vp J Q 

Finally, introducing z — cos0 [so that H(z) = (m/27r)(l — z)], we obtain 

SxV°(Q, 0) = f dz (5z 4 - 3z 2 ) (1 - z). (4.17) 

The relevance of the non-analyticity in the polarization bubble is now transparent: if Tl(z) was z-independent, the 
integral over z would vanish. However, because of the non-analyticity, II(z) varies linearly with z. The integral over 
z then does not vanish, and performing the integration we obtain 

«OHxoA(^) 3 M, (4. 18 ) 

where xo = 211(0, 0) = m/n is the static susceptibility of noninteracting fermions. 
Using (|4.10() . we then obtain the total contribution of diagram 1: 

_oW*=<V^_„ 4 f™U\ 2 \Q\ 



6 X i(Q,0)=25 X r(Q,0)= Xo - (4.19) 

Diagram 2 is another self-energy insertion into the particle-hole bubble. For a contact interaction, 5x2 is exactly 
(—1/2) of Sxi, the rescaling factor —1/2 comes from the fact that compared to diagram 1, diagram 2, has one lass 
fermionic loop with more than one vertex, and lacks the factor of two due to spin summation. Therefore 

2 /mC/ x2 



**,(«,()) = -*> j£ (4-20) 

The next diagram, diagram 3, represents a vertex correction to the particle-hole bubble. The q — contribution 

to this diagram can be shown to be of the same magnitude but opposite sign as the q = part of diagram 1 . To see 
this, we write the q = contribution to diagram 3 as 

SxT°(Q,0) = -^j HI d2fc ( 2 y^ Go(k,c,)G (k + q,c, + ») (4.21) 

xG (k + Q + q, uj + n)G Q {k + Q, w)n(?, O) (4.22) 

and consider a combination 

C=\s x r + Sxr°- (4-23) 
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Linearizing the fermionic spectra according to Eq. (|4.11|l , we re- write C as 



(4.24) 



where 



and 



Si = f de k G 2 {k, u;)G (k + Q, ^)G (k + q, lu + fi) (4.25) 



S 3 = J de k G {k 7 iv)G Q {k + q,uj + fl)G {k + Q + q, u + ^)G (k + Q, u). (4.26) 
Integrating over e k in Eas. H4.25l4.26f) yields 



Si = -27r£sgn (Q) (w (O - w)) - T 



S 3 = 2magn(n)9(u(Q-u})) 



in + v F k • q ) in + y p k ■ q-^ k ■ Q 
1 1 



vpk ■ Q iQ + vpk ■ q 
1 1 



JQ + vpk ■ q— vpk • Q iCl+ vpk ■ Q, 
Adding Si and S3 and performing some elementary transformations, we obtain 

Si + S 3 = 27risgn (CI) 6 (u> (Q - w)) ' ' 



(iSl + vpk ■ q) 2 ^ + l 'F k ■ V+v F k • Q 
Substituting the last expression back into Ea. (|4.24|l and making the change of variables k — > — k, q — > q results in 



C = - l -8 X q r\ (4-27) 



Together with (|4.23[1 . this proves that 



5xT°(Q,°) = Sxr°(Q,0) = -xo I 0. (4.28) 

The 2fci?-contribution from diagram 3 must be computed independently. The computations are performed along 
the same lines as for diagram 1. We present them in the Appendix El We obtain 

(«0 

Comparing this with Ea. H4.20[) . we see that, for a constant interaction, the 0(|<3|) contributions to diagram 3 from 
the singularities at q — and q = 2kp cancel each other. This result appears to be quite general (the same is true 
for D = 3 and D — 1 also (see below), but we do not know how to prove it other than to explicitly compute the 
diagrams. 

Next we consider diagram 4, which is obtained by inserting the particle-particle bubble into the original particle-hole 
bubble. Expressing 8\i via the product of four Greens's functions and the particle-particle bubble, we obtain 

MQ,o) = JJJ /^|^ Go (k,.) Go (k + Q,.) 

xG (q - k, fi - w)G (q -k - Q,fi- u)Tl pp {q, fi), (4.30) 

where H pp (q, fi) is given by (|2.8|) . 
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FIG. 4: A reduction of 5xi to the 2Jzf part of 5x3- 



In principle the result for 5 X a can be found by substituting the particle-particle propagator into (|4.H()|I . However, a 
straightforward approach is very cumbersome in this case. There is a more elegant way to compute 5 X i as the non- 
analytic part of this diagram is related to the non-analytic 2kp contribution from diagram 3, which we have already 
found. Indeed, it is easy to make sure that a non-analytic (oc |Q|) contribution from diagram 4 comes from internal 
momenta for which one of the internal 3-momentum transfers is small. We can then label the internal momenta in 
diagram 4 as shown in Fig. 21 an d set 3-momentum q to be small (there is a combinatorial factor of 2 associated 
with this choice). We can then represent diagram 3 as an integral-over-q of a product of two terms ("triads") each 
containing a product of three Green's functions: 

5x4 = -2x2U 2 J J -^-1 (q, O; Q) / (-q, -ft; -Q) , (4.31) 

where a "triad" is defined as 

I (q, O; Q) = / / — — 3-G (k, u) G (k - q, w - Q) G (k + Q, w ) . (4.32) 
J J (2tt) 

An extra overall factor of —2 in l|4.33|l is due to spin summation and the presence of one closed fermionic loop. At the 
same time, we can use the fact that in the 2kp part of diagram 3, one of the two momenta in the internal particle-hole 
bubble is close to incoming ones. Using the labeling as in Fig. 0] we can express the 2k p part of diagram 3 as 

5xf F =^J j^[/(q,tt;Q)] 2 , (4.33) 

Carrying out integrations over ek and ui in Ea. (|4.32|l . we find that 

J(-q,-n;Q)=-J(q,0;Q), (4.34) 

and hence 

Sx i (Q,0) = 6 X ^(Q,0) = X o^ {^j ^- (4.35) 

Similarly, diagram 5 differs by a factor of —1 from diagram 3 (the lack of the spin factor of two, compared to 
diagram 3, is compensated by an extra combinatorial factor of two). For a contact interaction, the non-analytic part 
of this diagram vanishes in the same way as it does for diagram 3. 

Finally, for the charge susceptibility, diagram 6 just differs by —1 from diagram 3, and diagram 7 differs by an 
extra —2 from diagram 4. For diagram 6, the extra —1 is due to the fact that, compared to diagram 3, q — and 
q = 2kp contributions are interchanged. For diagram 7, the extra factor is due to the spin summation and reflects 
the presence of two closed fermionic loops in diagram 7, as opposed to one loop in diagram 4. 
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Collecting all terms, we obtain 



Sx2(Q,o) 
Sx 5 (Q,o) 

As a result, 



§xi(Q,o) 
1 



X0 37 



mU 

4-7T 



<$xi(Q,o), <Sx3(Q,o) 



0, 6 X 4(.Q,0) 



1 



<*xi(Q,o), 



o, s X6 (Q,o) = o, s X 7(Q,o) = -Sxi(Q,o)- 



(4.36) 
(4.37) 



<5 Xc (g,o) = o 



mU 

4-7T 



(4.38) 



This result is consistent with the conjecture by BKV, who found that the spin susceptibility has a Q 2 In |Q|-dispersion 
in 3D, and conjectured that Xs(Q,fy should scale as |Q| in 2D. We emphasize, however, that we present for the 
first time an explicit calculation of Xs(Q,0) in 2D. BKV did not explicitly consider the charge susceptibility, but the 
absence of the non-analytic momentum dependence of Xc can be readily extracted from their analysis. 



2. D = 3 and D = l 



For completeness, we also performed full calculations in D = 3 and D = 1. In both cases, the results, Sxl D (Q, 0) cx 
Q 2 InQ, S\l D (Q, 0) oc InQ, have logarithmic non-analyticities in Q , which allows one to expand in Q from the very 
beginning. Doing so, we reproduced the results by BKV. 

In 3D, we obtained for the spin susceptibility 



5X3(Q,0) = 6 X5 (Q,0) = 0;<5 X2 (Q,0) = -^Xi(Q,0);5 X4 (<2,0) = ^ Xl (Q,0); 



5 Xl (Q,0) = 2S x r°(Q,0) = ±xl D (^ 



Q_ 



In 



Q 



(4.39) 



where Xo D = T^kp/TT 2 is the static spin susceptibility and a = mU/An is the scattering length. Combining all 
contributions we obtain 



kp 



In- 



(4.40) 



Eqs. H4.39fl an d l|4.40|l precisely coincide with the earlier results by BKV |2g. We also considered the charge 
susceptibility and found that, as in 2D, it does not possess a non-analytic dependence on Q. 
In ID, the relations between various components of Sxl D (Q,0) are the same as in 3D, and 



Sxl D (Q,0) = <5xi(Q,0) = 25 X r°(Q,0) = -2x1° C^j' 



k F 
m -^r- 



(4.41) 



This form 5x1 {Q,0) agrees with the well-known result by Dzyaloshinskii and Larkin |5£ 



B. Spin and charge susceptibilities at finite T and Q — 



In this Section, we consider the uniform (Q = 0) spin and charge susceptibilities at finite T . Of particular interest 
here is the question whether a non-analytic momentum dependence of the static susceptibility at T = is accompanied 
by that of the static susceptibility. We remind that in D = 3, according to Carneiro and Pethick |43 and BKV, 
x(Qj 0) — x (0, 0) behaves as Q 2 In |Q|, but x(0, T)~x (0, 0) is analytic and behaves as T 2 . Misawa |4l|, on the contrary, 
did find a T 2 InT-behavior. BKV conjectured that for a generic D, the momentum and temperature dependences of 
Xs should have the same exponents. 

As it was pointed out in the Introduction, there were two microscopic calculations of x(0,T) in 2D: by BKM 
[3fj| and CM .31]. Both groups found Xs(0,T) cx T and associated this non-analytical T dependence with the 
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square-root singularity in the quasiparticle interaction function / (k, k') caused by 2Uf scattering. We recall that 
the quasiparticle interaction function, / (k, k') , is obtained by computing the vertex T (k,w; k', u>'; q,0) to the second 
order in the interaction and using the relation pjj.[Tll] 

/ (k, k') = AT (k,e fe ; k'e k ,; q/Q -> 0) , (4.42) 

where A is a normalization factor, BKM [3(j explored the singularity / (k, k') at T — and for small but finite 
quasiparticle energies, efc and ek>. In their approach, the T— dependence comes from the Fermi functions. In the 
diagrammatic language, the approximation made by BKM accounts to evaluating the particle-hole polarization bubble 
near 2Uf at T — but at a finite frequency. CM included this effect into their consideration, but they also took into 
account a \T— singularity associated with the thermal smearing of the 2/c^-feature in the susceptibility. 

We compute Xs(Q = 0, T) in a straightforward diagrammatic approach (the same we employed for the case of 
Q ^ 0, T = 0), in which all possible sources of T-dependence are taken into account automatically. Our result differs 
by a factor of 2 compared to that of CM. We could not establish the reason for the discrepancy. 

We first report our results for D = 2 first and then analyze the case of arbitrary D. 



1. D = 2 



The analysis of x(0, T) proceeds in the same way as in Sec II V Al We found that the interplay between the non- 
analytic terms in various diagrams for the susceptibility at T ^ is exactly the same as at T = 0. Namely, the non- 
analytic pieces originate from the q — and 2kp non-analyticities in the particle-hole susceptibility, or alternatively, 
from the q — non-analyticity in the particle-particle susceptibility. We explicitly verified that the relative coefficients 
between non-analytic terms are the same as at T = 0. This implies that (i) just as at T — 0, there is no non-analytic 
T dependence in the charge susceptibility, and (ii) to obtain the full correction the spin susceptibility, it is sufficient 
to evaluate just one non-analytic contribution, e.g., <5x^ -0 (0, T). The full 5x s (0,T) is then given by 

6 X s(0,T) = 26xr°(0,T). (4.43) 

At finite T and Q = 0, a general form of <5xf -0 (0, T) is 

<5xr°(0,T)--8(7 2 r 2 ]T f f^^-Gl(k,u})Go^ + ci,oj n + fl m Mq,il). (4.44) 

Expanding the quasiparticle spectra near the Fermi surface, integrating over and then evaluating the sum over ui n , 
we obtain, after simple algebra, 



where Xq D = rn/ir, 



•i\T"m.ri: J.,?" (^) 2 I(T), (4.45) 



w-W-W*? 1 - ( " 6) 



and x = vpq. Expression (|4.4t)|) is rather tricky, because I(T) is formally ultraviolet-divergent. The most straight- 
forward way to get rid of the ultraviolet divergence is to introduce a short-range (lattice) cutoff in the momentum 
integral so that x < Xq ~ W. Evaluating the integral over x first we obtain 



where 



S (m) = 1 + 2 ""V, - 3 ^ 5 . (4.48) 
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For Vl m <C X , S(m) is close to 1, i.e., S (m) = 1 + 0((r^/ 'A ) 2 ), whereas for Q m > X it falls off rapidly [as 
(X /n, n ) 2 }. The vanishing of S(m) at large m ensures the convergence of the sum in Ea. ()4.47|) . and allows one to 
use Euler-Maclaurin formula Applying it, the sum reduces to 



T 

AE F 



OO rp 

n— — oo 



dmS*(m) — S*'(0) + 

v ; 2AE F K 1 



2E F 



(4.49) 



where . . . stands for higher-order derivatives of S. All derivatives of S(m) obviously vanish in the continuum limit 
W — > oo. The remaining integral term in (|4.49|) gives 



T 
2E~ f 



dmS (m) = — . 

o K ' 16 E F ' 



(4.50) 



which is a T- independent contribution. As a result, the above computation does not yield a linear-in-T piece in 
<5xr>,T). 

A more careful inspection of the steps we took to arrive at this result reveals a problem. Namely, it is obvious 
from H4.46fl that the term with m — 0, i.e., with Q m — 0, vanishes for any finite q. However, in the sum in (|4.47|) 
the m = term is present and contributes T/AE F . As the static susceptibility is properly defined as the limit of 
x{Q,T) at Q — > 0, one should always keep q finite at the intermediate steps of the computations. Alternatively, one 
can perform calculations for a finite system and then extend the system size to infinity. In both cases, there exists a 
lower cutoff in the integral over q. This cutoff plays no role for all terms with m ^ but it eliminates the term with 
m = 0. Subtracting off this term from ()4.47|) . and using our previous results we obtain a universal, linear-in-T piece 
in J(T) 



J(T) = - 



T 
4E^' 



(4.51) 



An alternative way to arrive at Ea. (|4.51J) is to perform the summation over fl m in l|4.46(l first, keeping q finite, and 
then integrate over q. Performing the summation, we obtain 



J(T) 



1 



4E F J 



dy y 



' dy 1 



1/2 



l/2\ 



I*- 

dy 2 



3/2 



riBiy 



1/2N 



(4.52) 



where tib (z) = (cxp (z/T) — 1) 1 is the Bose distribution function and y = v 2 F q 2 . Integrating by parts, we obtain 
from (|4.52(l 



I(T) 



■—(! + - 
2E F V 2T 



' dy^-{y^n B {y^ 2 )} 



T 
1~Ef~'- 



(4.53) 



in agreement with l|4.51|l . 

The above analysis shows that Xs(0,T) does indeed contain a linear-in-T term in D — 2, as earlier studies conjec- 
tured. However, the physics behind this term is very different from the one that leads to the \Q\ piece in x s (Q,0). 
Substituting I4.51|) into I4.45[> and then using l|4.43l) , we obtain 



5 Xs (0,T) = 2 X t» I — ) — . 



mil 



An 



T 



(4.54) 



This is the central result of this subsection. 

We remind the reader that the full x s (0,T), given by (|4.54(l comes from the dynamical particle-hole bubble. 



To 



emphasize this point, in Appendix [0 we compute 5xi F neglecting the frequency dependence of the polarization 
bubble, and show that this yields an incorrect prefactor in the linear-in-T piece. 

We did not attempt to verify our Sxi by explicitly computing a linear in T contribution from 2k F polarization 
bubble at a finite T (as we did for \Q\ term at T = 0). This calculation would require, as an input, the analytical 
expression for the dynamical polarization bubble near 2k F at a finite T. We couldn't obtain this expression in 
a manageable form, nor we could find it in the literature. It would be interesting, however, to verify our Sx1 kF 
numerically by using the numerical results for H(q,u>,T) |6lj |. 
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2. other dimensions 

For arbitrary D, the consideration analogous to the one for D = 2 yields, instead of Eq. I]4.45[l . 

5 X d(0,T) = -CU 2 I D (T), (4.55) 

where C is a positive constant, 

I D (T) = l^T^ / a***- 1 (n2 % )3 ( 4 - 56 ) 

rn \ m J 

and x = v F q. For D = 2, jO^I coincides with (|4~4ll modulo a piece [T^ m /fa (^-^ 2 ) /(^m + ^ 2 ) 3 ] that 
vanishes upon integration over x. The ambiguity with the order of summation and integration was resolved in the 
previous section; now we know that it is safe to sum over Q m first and the integrate over x. Performing the summation 
with the help of the well-known formula 

1 l + 2n B (x) 

^ ^LTx~- = 2x ' (4 ' 57) 



we find 



/d(T) = 2E F T=l ' dXX 



D-l 



l + 2n B (x) i ^ 2 d l + 2n B (x) | x 4 d 2 l + 2n B (x) 



2x dx 2 2x 2 dx 2 2x 



(4.58) 



Evaluating the integral over x and introducing an infinitesimally-small 5 to eliminate infrared divergences at interme- 
diate steps, we find the T— independent part of Id (T) for D > 2 to be given by 



D—l ^oo i n-2 



8 \E F J J e z - 1 

= J D - 2 ^ A - D ) T (D-l)C(D-l), (4.60) 

where T (x) and C (x) are the Euler and Riemann functions, respectively. For D — » 2, the pole of the £ -function, 
C(Z) — 1) — * 1/(D — 2), is canceled by the prefactor D — 2, so that Ii^T) is finite and equal to —T/AEp, in agreement 
with lj433j) . 

For D < 2, care has to be taken to ensure the cancelation of the divergent terms. The final result for this case is 

(2 — D)(A — D) ( T \ D1 f°° dz (I 



im-^g — ^ i («') 

We see that for arbitrary D, the function Id(T) (and thus the spin susceptibility) scales as T D ~ 2 . In an explicit form, 



<$X(0,T) = - 




(4.62) 



Function / (D) diverges logarithmically for D — l (and at D = 1, S\ oc InT). Near D = 3 function f(D) is perfectly 
regular and equal to 

/(3) = ~ (4.63) 

As we see from Eas. (l4.62|l and (|4.63() . this last result implies that in 3D, the leading temperature correction to 
the susceptibility scales as T , and there is no logarithmic prefactor. This agrees with the results of Carneiro and 
Pethick and BKV. 

Obviously, the absence of the T 2 In T-behavior of x(0,T) in 3D, and Q 2 In Q-behavior of x(Q, 0) implies that there 
is no one-to-one correspondence between thermal corrections and quantum corrections at finite T. Our consideration 
indeed shows that thermal and quantum corrections are not equivalent. 

We also see that although f(D) goes smoothly through D = 2, the functional form of f(D) changes between D > 2 
and D < 2. The consequences of this fact are, however, unclear to us. 
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k k+q -k+1 -k+m 



FIG. 5: One of the self-energy diagrams for (7(0) 7^ U(2kp). Momenta q, I and m are small compared to k. 



V. FINITE-RANGE INTERACTION 



In the previous sections we considered the model case of a contact interaction, characterized by a single coupling 
constant U which is independent of the momentum transfer. Now we analyze the more realistic case of a finite-range 
interaction when the coupling is a function of the momentum transfer U — ► U (q) , where U (q) is such that U(0) and 
U(2k F ) are finite. Our key result is that only these two parameters are important. 



A. Self-energy 

We begin with the self-energy. For momentum-dependent U(q) the two self-energy diagrams in Fig. n uave to be 
considered separately. For diagram shown in Fig. [IJ,, the extension to U = U(q) is straightforward-the factor 2U 2 
for that part of the self-energy which corresponds to process b) in Fig. 2 (we recall that only that part contributes 
to thermodynamics) is be replaced by U 2 (0) + U 2 (2k F ). The diagram in Fig. ^> requires more care, but we know 
from the analysis of the "sunrise" diagram for the self-energy (Fig. 4b) that a non-analytic piece comes from the range 
where two internal momenta in the self-energy diagram are near — k, and the third is near k. For diagram Fig. ^d, 
this implies that the momenta are labeled as in Fig. El 

It is then obvious that the overall factor for the diagram in Fig. QJd is U(Q)U(2k F ). Process a) in Fig. 2 determines 
that part of the self-energy which is singular on the mass-shell and does not contribute to thermodynamics. The 
overall factor for that part is U (0) 2 . Collecting all contributions, we find that 



mU 2 (Q) , 2 W 
w > T = 1fi 3 2 w + 7rT ) ln I 1 



T2 



m(U 2 (0) + U 2 (2k F )~U(0)U(2k F )) ( 2 2 A i . ( u \\ ( , u 

where A and A are constants, and the scaling function f(x) is given by (|3.25l) . The real part of the self-energy is 
given by 

m([/ 2 (0) + [/ 2 (2fc F )-t/(Q)t/(2fc F )) , U \ 

= ^2 "H« (t) ■ ( 5 - 2 ) 

The limiting forms of the scaling function g{x) are 5(00) = 1 and g{x <l)w Alog2 / (x) . 



B. Spin and charge susceptibilities 

The same consideration holds for the susceptibility-the very fact that all non-analytic contributions come from the 
vertices with near zero total momentum and transferred momentum either near zero or near 2kp implies that for 
U = U(q), an overall factor of U 2 is replaced either by U 2 (0) or U 2 (2kp), as in diagrams 1, 3, 6 and 7 in FigEI an d 
by U(0)U(2kp), as in diagrams 2, 4 and 5. With this substitution, we have, finally 

6 X i(Q,T) = K(Q,T)(U 2 (0) + U 2 (2k F ));S X 2(Q,T) = ~K(Q 1 T)U(0)U(2k F ) 

5 X3 (Q,T) = K(Q,T)(U 2 (2k F )-U 2 (0));SXi(Q,T)=K(Q,T)U(0)U(2k F ); 

6 X5 (Q-T) = 0, S Xe (Q,T) = ~K(Q,T)(U 2 (0) - U 2 {2k F ))- 

SX7(Q-T) = -K(Q,T)(U 2 (0) + U 2 (2k F )) 7 (5.3) 
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FIG. 6: Another way of labeling momenta in 5x5- 



where K(Q,0) and K(0,T) are given by Eqs. I|4.18|l and (|4.54|) . respectively: 



21) 2 ( mU\ 2 \Q\_ T s fnm ^_„.2D ( mU\ 2 T 



where xo = m/%. When both Q and T are nonzero, K(Q,T) — K(Q,0)g(vpQ/T), where g(x) is a scaling function 
subject to g(x 3> 1) oc 1/x. However, we did not attempt to compute g(x) at x other than x = and x = oo. 
Collecting all contributions we find for the spin susceptibility 

5 Xs (Q, T) = 2K(Q, T)U 2 {2k F ). (5.5) 

As for the case U =const, the charge susceptibility is regular because all non-analytic corrections from individual 
diagrams cancel out. 

Eqs. Q5.1|) . I|5.2|) . I|5.4|l and l|5.5|l are the central results of this paper. 

While it is intuitively obvious that the momentum dependence of the susceptibility should only include U (0) and 
U(2kp), this intuition is based on the analysis of the self-energy but not the susceptibility itself. It is therefore 
worthwhile to demonstrate explicitly that non-analytic terms in the susceptibility do not depend on the momentum- 
averaged interaction. This is what we are going to do in the remainder of this Section. 

To demonstrate that only U(2kp) matters, consider one of the diagrams for which, as we claim, the non-analytic 
term scales as U(0)U(2kp), i.e., diagrams 2, 4 and 5. Each of these diagrams has two interaction lines. Quite obviously, 
one of momentum transfers should be near zero. The issue is to prove that the other one is near 2k F . Consider for 
definiteness diagram 5. The net result for this diagram is zero, but this is a result of the cancelation between two 
contributions, 5\t{Q) an d 3xt(Q)> which differ in the choice of which of the two interactions carry small momentum. 
Consider one of the choices. We label the internal momenta in the diagram as k, k + q, k + Q, k + q + Q, 1 + q/2, 
1 — q/2, where Q is the external momentum, and introduce two angles 8\ and 82 between q and 1 and between q and 
k, respectively (cf. Fig©. 

The integration over k and the corresponding frequency u> is straightforward (see Appendix . Introducing then 
q = r cos <p and O = r sin </>, where f2 is the frequency associated with q, we integrate over r and, after redefinition of 
the variables, obtain that the non-analytic, linear-in-Q piece of diagram 5 reduces to 

where 

J = f°° dxx 2 r ddl - r 7 U ( 2k F sin 2 6 -^.) . (5.7) 

Jo J-n X + ICOS0! J_, {x + icos9 2 ) 4 V 2 J v ; 

For a constant interaction U(q) — U, we can integrate independently over 9\ and 62, and then integrate over x, 
which gives J — ir 2 /6. The result for 8\t{Q) then coincides with one of the two contributions to SxsiQ), as we 
discussed in Sec lIV ATI A relevant point here is that typical cos 6^2 are of order x, whereas typical x are of order 
1. Hence 81—82 ~ 1, i.e., typical angles between two momenta are generic. This implies that the argument of 
U{2kp sin 2 (6*i — 8 2 )/2) is just of the order of kp but not necessarily close to 2kp. 
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We now show that, in fact, only B\ — 82 = ±7r matter. To see this, we introduce diagonal variables a = (61 + 62) /2 
and b = (61 — 8 2 )/2 and integrate first over x and then over a. This integration is tedious but straightforward, and 
carrying it out we obtain, after some algebra, 



L 



7I-/2 

dbU(2k F sin 2 6)Re [S(b) + S(ir - b)} , (5. 



where 



Then 



S(b)= (^ + ^^)lncos26(— =■ : — -). (5.9) 

' 1 3 sin 4 b) V sin 2b - iS sin 2b + iS 



J = i51mf rizlncoszsin 2 z+S 2 (- + - C ° S f ) U{2k F {sm z/2) 2 ) (5.10) 
J a V 3 (smz/2) 4 / 

The integral does not vanish due to divergences near z = and z = ir. The divergence near z = does not contribute 
to the imaginary part of the integral, but the one near z = n does contribute. Restricting z near ir, we obtain 

j= i u{2kF) r^=Y u{2kF) - (5 - n) 

This consideration shows that although for a momentum-independent interaction we could evaluate Sxt(Q) in a scheme 
in which the angular integrals were not restricted to a particular 8\ or 62, the calculation performed in another way 
demonstrates that the whole integral comes only from the range where Q\ — 8 2 = ±7r. For a momentum-dependent 
interaction, this implies that only U(2kp) matters, precisely as we anticipated. Similar calculations can be repeated 
for other cross diagrams with the result that the overall factor is always {7(0){7(2/cp)- 

The above consideration is another indication that the non-analyticities in the specific heat and spin susceptibility 
come from the two interaction vertices in which the transferred momentum is either near or 2kp, and simultaneously 
the total momentum for both vertices is near 2kp. 



VI. CONCLUSIONS 



We now summarize the key results of the paper. We considered the universal corrections to the Fermi-liquid forms of 
the effective mass, specific heat, and spin and charge susceptibilities of the 2D Fermi liquid. We assumed that the Born 
approximation is valid, i.e., mU{q)/4Tr <C 1, and performed calculations to second order in the interaction potential 
U{q). We found that the corrections to the mass and specific heat are non- analytic and linear in T, and obtained 
for the first time the explicit results for these corrections. We next found that the corrections to the static spin 
susceptibility are also non-analytic and yield the |Q|-dependence of Xs(Q, T — 0) and T dependence of Xs(Q = 0, T). 
We obtained for the first time the explicit expressions for the linear-in-Q and linear-in-T terms in the susceptibility. 
We found that the corrections to the charge susceptibility are all analytic. We also performed calculations in 3D and 
confirmed the results of BKV and others that the correction to Xs{Q, T — 0) scales as Q 2 InQ, but the correction to 
Xs(Q = 0, T) scales as T 2 without a logarithmic prefactor. 

We analyzed in detail the physical origin of the non-analytic corrections to the Fermi liquid and clarified the 
discrepancy between earlier papers. We argued that the non-analyticities in the fermionic self-energy and in Xs{Q, T) 
are due to the non-analyticities in the dynamical particle-hole susceptibility. We argued that the non-analyticities 
in the fermionic self-energy and in Xs(Q,ty ar e due to the non-analyticities in the dynamical two-particle response 
functions. We have shown that non-analytic terms in the self-energy and the spin susceptibility come from the 
processes which involve the scattering amplitude with a small momentum transfer and a small total momentum. We 
explicitly demonstrated that the non-analytic terms can be viewed equivalently as coming from either of the two 
non-analyticities in the dynamical particle-hole bubble-the one near q = and the other one near q = 2kp-or from 
the singularity in the dynamical particle-particle bubble near zero total momentum. We also demonstrated explicitly 
that the non-analytic terms in all diagrams for the susceptibility and the self-energy depend only on U(0) and U {2k f), 
but not on averaged interactions over the Fermi surface. Only under this condition, is there a substantial cancelation 
between different diagrams for the susceptibility. Due to these cancelations, the non-analytic correction to the spin 
susceptibility depends only on U{2kp), but not on {7(0), and scales as {7 2 (2fcp)- The non-analytic corrections to the 
effective mass and the spin susceptibility scale as {7 2 (0) + U 2 {2kp) — U{Q)U{2kp)- 
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-k+q - k -q' 

b) 

FIG. 7: Two processes contributing to the log-singularities in the self-energy. 

The non-analytic Q behavior of x s (Q,T = 0) obtained in both 3D and 2D questioned the validity of the Hertz- 
Millis-Moriya phenomenological theory of quantum phase transitions. This theory assumes a regular ^-expansion of 
the spin susceptibility. Indeed, extending the results for Xs(Q, T) to the critical region one obtains a rather complex 
quantum critical behavior 62], which is very different from the Hertz-Millis-Moriya theory. We caution, however, 
that the non-analytic behavior of \s{Q,T) was obtained within the Born approximation, when fermions behave 
as sharp quasiparticles. Near a magnetic transition, the fermionic self-energy is large, and destroys the coherent 
Fermi- liquid behavior beginning at a frequency which vanishes at the quantum critical point. In this situation, the 
second-order perturbation theory is unreliable. The issue whether non-analytic corrections to the static Xs{Q,T) 
survive at criticality is now under consideration and we refrain from further speculations on this matter. 
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APPENDIX A: MASS-SHELL SINGULARITY 

In this Appendix, we take a deeper look into the origin of the logarithmic divergence of the self-energy on the 
mass shell. To better understand where it comes from, we come back to the derivation of (|3.11|) . Re- writing (|3.11|) 
as (I3.12all to logarithmic accuracy, we now argue that the two logarithmic terms in i|3.12a[l come from two different 
processes, as shown in Fig[7| In the first process (FigEJi), all four momenta are close to each other, and in the second 
one (Fig[7p), the net momentum of the two incoming particles is close to zero, whereas the momenta of the outgoing 
particles are close to their initial values. In terms of the momentum transfers, both processes are of forward-scattering 
type. To see this, we notice that for generic u)/ek, i.e., not too close to the mass shell, the logarithmic form of the 
self-energy is due to 1/Q behavior of the momentum integrand at vfQ 3> £l,tv. This l/Q form in 2D results from 
the combination of two facts: (i) the polarization operator n(Q,f2) behaves as tt/vpQ, and (ii) the imaginary part 
of the fermionic propagator, integrated over the angle between Q and external momentum k, behaves as l/Q. The 
product of the two terms yields J QdQ/Q 2 that gives rise to a logarithm. It is easy to make sure that for VfQ S> fi, 
typical values of are close to ±7r/2, the deviation from these values being of order |f2| /vpQ. That means that the 
external momentum (k) and the internal (small) one (Q) (as labeled in Fig^Jt) are nearly orthogonal to each other. 
The same reasoning also works for the polarization bubble. If the two internal momenta in that bubble are p and 
p + Q, then typical p and Q are also nearly orthogonal. Since both k and p are orthogonal to the same Q, and both 
are confined to the near vicinity of the Fermi surface, they are either near each other, or near the opposite points of 
the Fermi surface. If p and k are close to each other, all three internal fermionic momenta in the second-order diagram 
are close to external k, if p is close to — k, out of three internal momenta one is close to k, while the other two are 
close to — k. These two regions of intermediate momenta give rise to two logarithms in i|3.12al) . The logarithm that 
diverges on the mass shell comes from a region where all momenta are close to k. To see this, we recall that the actual 
divergence is the consequence of the fact that both the polarization bubble and the angle-averaged G"(k + Q, lu + ft) 
at the mass shell possess square-root singularities in the form 1 / y/ (vfQ) 2 — (f2) 2 such that the product of the two 
gives (vfQ) 2 — (^) 2 ) ; and the momentum integral diverges. The square-root singularities come from near parallel 
p and Q and k and Q, respectively. Obviously then, k and p are near parallel, i.e., they are located near the same 
point at the Fermi surface. With a little more effort, one can show that as u> approaches e&, typical angles between 
p and Q and between k and Q, both move from near 7r/2 (or — 7r/2) to near zero, but in such a way that k and p 
remain parallel. This once again confirms that the divergent logarithm comes from the process in FigE^ ( a U internal 
momenta are close to k), while the "conventional", non-divergent w 2 lno;-term comes from the process in Fig|7p. 
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FIG. 8: Non-trivial diagrams for the self-energy in ID. ± denotes the propagator of a right/left moving fermion. 

The analysis can be extended to finite T, and the (anticipated) result is that £'/ given by l|3.22(l comes from forward 
scattering, while given by (|3.23f) comes from back scattering. 

It is interesting follow the same arguments for D = 1. In this case, processes in Figd acquire even simpler physical 
meanings: process a) is forward scattering of fermions of the same chirality, e.g., two right-moving fermions scatter into 
two right-moving ones, whereas process b) is forward scattering of fermions of opposite chirality, e.g., a right-moving 
fermion scatters at a left-mo ving one so that their respective chiralities are conserved. In the g-ology notations, vertex 
a) is (72 and vertex b) is gn [5l| . In the Luttinger model, when only forward scattering is taken into account, the 
self-energy of, e.g., right-moving fermions is represented by the set of diagrams shown in FigEl J52], where ± denotes 
propagators of right/left moving species 

G± (k, lo) = — 4 =v F (k T k F ) . (Al) 
iuj-e k 

Diagrams a) and c) contain two vertices of type a) from Figd whereas diagram b) contain vertices of type b) from 
Figd The imaginary parts of the retarded polarization bubbles for right- and left- moving fermions for \Q\ — > are 
given by 

riR ± =®5(Sl T v F Q). (A2) 

The delta- function form of H R ± is due to the fact that in ID and for |Q| — > the particle- hole continuum shrinks to 
two lines in the (fi, Q) plane described by O = v F \Q\ . The combination of the diagrams a) and c) in Fig[S|yields for 
the imaginary part of the self-energy 

[ S 'A + (k,-)] a)+c) = ^^(--4), (A3) 

F 

We see that given by l|A3jl . which is a ID analog of our EJ( from H3.12c|) . is very singular on the mass shell but 
vanishes outside the mass shell. At the same time, diagram b) in FiglHI yields 



^-(H- |4|), for m>|4| 
0, otherwise. 



This self-energy vanishes on the mass shell, but for a generic co/ek it yields (k, cj)] ^ cx |w|. This \lo\ dependence 

obviously implies that Fermi-liquid behavior is in danger. 

Which of the two terms is actually relevant? In ID, the answer is well known: the summation of infinite series of the 
diagrams yields the non-Fermi-liquid behavior, and the resulting state-the Luttinger liquid-is free of singularities on 
the mass shell. This implies that the mass shell singularity of Eq. i|A3|) is completely eliminated by the re-summation 
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of diagrams to all orders in the interaction. This can be shown explicitly either via Ward identities or using the 
bosonization j5^|. Furthermore, the exact solution of the model with only type a) scattering ("174 -model") yields a 
/ree-Fermi-gas behavior with a renormalized Fermi velocity, i.e., no mass-shell singularity. This all implies that the 
mass shell singularity found in the second-order self-energy diagram in ID is an artificial one and is eliminated by 
higher order diagrams. 

The same elimination of the mass shell singularity holds in 2D, as we now demonstrate. Indeed, as we mentioned 
before, the logarithmic divergence in l|3.12a|) at w — ek is the consequence of the matching of the two square-root 
singularities: one resulting from the angular integration of the fermionic Green's function, and another one being 
the 1/^(vfQ) 2 — ft 2 singularity in W R (Q,ft). Suppose now that the interaction gets renormalized (screened) by 
higher-order terms in U so that U -> U(Q, ft). The combination U 2 T1 R (Q, ft) in (£3 is now replaced by U&(Q, ft). 
In the RPA approximation (which is not a controllable one for a short-range interaction), 

U" (Q ft) = ggiM) 

[l + t^(Q,tt)] 2 + [C/ri£(tt,Q)] 2 



= 2^2 n V(vFQ) 2 -ft 2 0(v F Q-\ft\) 
m (1 + jjy [(v F Q) 2 - ft 2 } + U 2 ft 2 ' 

where U = mU/2n. Obviously, U R now vanishes at Q = \ fl\ /uf, and the divergence is eliminated. At the same time, 
the logarithmic dependence on ft in (|3.9|l . and hence the w 2 lnu; form of the self-energy, survive as they come from 
typical ft ~ vfQ for which H' R (Q, ft) and H R (Q, ft) are of the same order, and hence the screened interaction is of 
the order of the bare one. Note that this reasoning is also valid for the Coulomb interaction, for which the RPA 
approximation is asymptotically exact in the high-density limit. 

Another argument that the mass-shell singularity is artificial is that it is eliminated, already at the second order 
of interaction, if one takes into account the curvature of the fermionic dispersion. Indeed, in obtaining (|3- 1 XI) . we 
linearized the fermionic dispersion near the Fermi surface, i.e., approximated tk+q by €k + v F qcos9. Using the full 
quadratic dispersion, we obtain, instead of (|3.10|l 

ml] 2 W 2 
S«(k,w) = ■5-5-5- / dflftln ][2fi-o; + e fc ]+A(w,n)| , (A5) 

oir a v F J Q € k - U) 

where 

A (w, ft) = ^- (3w - e k - ft) (A6) 

and where, for the sake of definiteness, we assumed uj > 0. On the mass shell, u) = ek, the integration over ft yields a 
finite result 

Y!' R (k,u)\ u=ek = — ^^ln— . A7 
lbn^vp \w\ 

The crossover between Eas. H3.11|) and (| AT|) occurs when, inside the log in Ea. (|3.11fl . A (ui, ft) becomes comparable to 
the other term in the denominator, i.e., when 

\oj-e k \ ~u?jW. (A8) 

For \u> — £k\ ^> oj 2 /W, the leading asymptotic behavior of S'^(k, oj) is given by l|3.11|l and for \u> — efe| uo 2 /W it is 
given by (|A7(I . A general formula which interpolates between the two limiting cases might, in principle, be obtained 
but wc do not dwell on this here. Notice that T, R (k 7 uj) on the mass shell is by a factor of 3/2 bigger than its value 
on the Fermi surface, which means that, for fixed u>, the slope of S^(k, uj) as function of becomes steeper as the 
mass shell is approached. 

The same result can be also obtained by calculating the quasiparticle lifetime for T — which, by definition, is 
taken directly on the mass shell. For D = 2, the Fermi Golden Rule gives 

U 2 m r u r° rW/v F f r 

1 / T ( UJ ) = ^3-J dn J dtJ ' J dQQ J d6 J d0'S{fl-e k + e k - Q )S{fl-e p+Q +e p ), (A9) 



where e = ,u>' = e p , and 9,9' are the angles between k and q and p and q, respectively, and W is the ultraviolet 
energy cutoff. For linearized dispersion the arguments of the first and second delta-functions in (|A9|) reduce to 
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+ vpQ cos 9, n — vfQ cos 9', respectively. Each of the angular integrations yields a factor of 2/ \J (vpQ) 2 — and 
the integral over Q 



A 



W/v F 



dQQ 



1 



n\/v F 



{v F QY - W 



(A10) 



diverges logarithmically at the lower limit. To regularize the singularity, one must to keep the higher-order terms in 
£k-Q and £ p +q. On the mass shell, 



£k — £k-Q 

£ P - £p+Q 



Vp I 



2E F ) 



cos 9 



2m 



cos 9' 



91 
2m 



Now the integral over Q takes the form 



where 



dQQ 

5 
6' 



(v F Q) 2 - & - 5 J {v F Qf - n 2 



oj Q 2 



hp m 



(Alia) 
(Allb) 

(A12) 

(Al3a) 
(A13b) 



The lower limit in the integral is such that the arguments of the square roots are positive. The momentum integral 
is controlled by Q ~ |f2|/i>F. To logarithmic accuracy, one can then just replace Q by |fi|/f.F in HA13aJA13b|l . After 
this replacement, the momentum integration can be easily performed and gives 



A = 



In- 



E 2 F W 



2v% n 2 (ui + w') + n 3 ■ 



(A14) 



We next have to perform the frequency integration. It is easy to verify that, in the two integrals over frequency, the 
dominant contributions come from the regions f2 ~ u>' ~ u>. To logarithmic accuracy, one can then simplify A to 



3 ln W 
24 D ^ 



(A15) 



We also used the fact that Ep ~ W. Substituting this into i|A9l and performing frequency integrations we obtain 
finally 



3U 2 m 



t (lu) 87t 3 i;|, 



> 2 In 



W 



(A16) 



We see that 1/t(u) is finite-the only memory left about the mass-shell singularity for the linearized spectrum is the 
enhanced numerical prefactor. Identifying 1/t with 2E'^, we see that the results for 1/r and £"(a> = e^) coincide, as 
indeed they should. 



APPENDIX B: POLARIZATION BUBBLE NEAR 2k F 



In this Appendix, we show that the computation of a non- analytic piece in the particle- hole bubble at Q « 2k F 
can be always performed in such a way that the dominant contribution comes from fermions near the Fermi surface 
and with nearly antiparallel momenta ±Q/2. We do this in two ways. First, we compute H p h{Q,Q m ) explicitly 
and check where the non-analyticity comes from. Second, we compute H p h(Q,Q m ) by linearizing the dispersion 
of fermions, forming the polarization bubble, near ±Q/2 and show that the non-analyticity in TlphiQ, Q m ) comes 
from the lower limit of momentum integration and therefore does not depend on the upper cutoff imposed by the 
linearization procedure. 



32 



1. Explicit computation 



Consider first T = 0. Labeling the momenta of internal fermionic lines in the polarization bubble as p ± Q/2 and 
for T — 0, we obtain in Matsubara frequencies 



n(Q,o r 



d 2 pduj Q _ . Q . 

(27r) 3 ^ P ~2~' ^ G (P _ Y'^n)- 



For a circular Fermi surface 



e p±Q/2 



P 



kp _j_ pQ cos | Q 2 



2m 2m 8m 

Substituting (|B2(I into IjBljl and integrating over frequency and then over p, we obtain for Q < 2kp 



mo n ) = — 1 1 - 2^ 

UW,iW 27T I 7rQ 2 



^ 2 dfl 



cos z 



arctan 



Pi 



arctan 



P2 

mf2 r 



(Bl) 



(B2) 



(B3) 



where 



Pi, 2 = Qcost 



kp 



— sin 2 d±-Q 2 cos 2 (9. 
4 2 



For Q = 2kp, we have pi = 4fcj, cos 0, P2 = 0, and l|B3(l reduces to 



n(Q,o m ) = 



2?r 



1 - 2- 



ttQ 2 



^ 1 2 d0 



■ arctan 



cos z 



4fc| cos 2 ( 
mO m 



m 
2^ 



1 - 



i /i a. 



l/2\ 



(B4) 



(B5) 



It is easy to see that the integral comes from cos 2 8 ~ \il\/Ep, i.e. typical p are nearly orthogonal to Q. Furthermore, 
in the integral over p, typical p were of order Qcos8. Hence typical p are of order y^rn|?l^J, i.e. at vanishing tt, the 
integration is indeed confined to internal momenta which nearly coincide with ±Q/2. 

The same reasoning is valid also for Q in a narrow range near 2kp. For Q < 2kp, Eq. 1B5|) can be re-written as 



n(Q,n m ) 



2- 



^/ 2 d# 



arctan ■ 



Q 2 cos 2 



cos z 



(mfi m )(l-Q 2 cos 2 8e 2 ) 



(B6) 



where e 2 = (Q 2 /4 — kp) /m\Sl m \. Assuming that the integral is dominated by 8 near 7r/2 and expanding 8 to linear 
order near 7r/2, we obtain after simple manipulations 



27T V irkp 



dz arctan - 



(B7) 



We see that the integral is convergent, i.e., the linearization of cos8 near ir/2 does not lead to cutoff-dependent 
integrals. This implies that the non-analytic piece in the polarization operator comes from typically small cos 8 and 
hence from typically small internal p oc cos#. Evaluating the integral over z in l|B7|) . we obtain 



777 

n(Q,fi m ) = — 



i 

i - - 

2 



V 



10, 



E F 



1/2 




(B8) 



where Q = Q — 2kp. This is the result that we cited in the text (Eq. J22})- 

For Q > 2fc_F, i.e., Q > 0, the calculations proceed in the same way. Integrating over p and over u> and again 
expanding to linear order near 8 = n/2 we obtain after straightforward manipulations 



n(Q,n m ) 



2tt 



Q_ 

kp 



1/2 



■Ky/2 V E F 



10, 



1/2 f l/2e 



dz arctan 



(l-4e 2 z 2 ) 1 / 2 
z 2 + e 2 



(B9) 



Evaluating the integral we find that the result reduces to Eq. (|B8 
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2. Another way of calculating H(Q ~ 2if,fi m ) 

For completeness, we also compute the non-analytic part in Tl(Q, f2 m ) near 2kp by explicitly restricting the integral 
over p in (|B1|I to small p and assuming that p is nearly orthogonal to Q. This calculation shows in a more direct 
way that typical values p are indeed small. To avoid lengthy calculations, we assume that Q = 2kp and aim at 
reproducing the yOm non-analyticity. For Q = 2kp, the energies on the internal fermionic lines are Ck F n+p and 
t-kpn+p- Introducing x — vpp and 7 = l/(2mv F ) — l/(4Ep), expanding cos# « 9, where 9 = it/2 — 9 and 
substituting into IjBlfl . we obtain 



1 

4ir 3 v F J_ oc 7_oo "~ Jo " (x6 + jx 2 - i(J n ) (x9 - jx 2 + i(uj n + fi m )) ' 



n(2fc F ,O m ) = ^-g-y / dfl / dw / — . — ^ — 77- (BIO) 



Introducing y = x9 and integrating over y, wc obtain after simple manipulations with variables 

1 f°° , f°° zdz 



n(2^,o ro) = ^y o d xj ™. (bid 



The integration is elementary and yields 



-1 /■ CXD j 

n(2fc F ,fi) = -^— / 4=. (B12) 



The divergence of the integral at the upper limit simply reflects that a constant term in the polarization bubble 
cannot be reproduced this way. However, the lower limit of the integral over z yields a universal and non-analytic 
contribution to H(2kF, 0, m ) of the form 

1 /|OJ\ 1/2 m/|ft' xl/2 



n smg(2fcF ,a„H-^^j =-to[i£) ■ (B13) 

This coincides with Eq. (|B8|I . It is essential that this result does not depend on the upper limit, and hence typical 
internal momenta scale with external fi. This obviously implies that typical values of p are indeed small. 



3. Finite temperature 



At finite T, a sharp ^JQ — 2kF non-analyticity in the static polarization operator is softened in qualitatively the 
same way as it is softened by a finite f£ m at T = 0. In general, 

m ( ( T \ 1/2 fv F \Q-2k F \ Vt m \\ 

n(Q,n^) = -^i-(-) rp F \ ^))- (Bi4) 

We could not find a simple analytical expression for the scaling function ^(a;, y) at arbitrary values of its arguments. 
At Q = 2k F and fi = 0, $(0, 0) w 0.339. 



APPENDIX C: EQUIVALENCE OF Q = AND Q = 2A" F CONTRIBUTIONS TO THE SELF-ENERGY 

In this appendix, we explicitly compute the contribution to the self-energy from the 2kp non-analyticity in the 
particle-hole bubble, and show that it is equal to the backscattering part of the self-energy from the q = non- 
analyticity. We will also show that the non-analytic self-energy can be equally viewed as coming from the singularity 
in the particle-particle channel at zero total momentum and frequency. 



1. 2&f part of the self-energy from the particle-hole channel 

Since our goal is to verify a general reasoning that q = and 2kp contributions to Y,(k,ui) are equal, we focus on 
the case T = e& = 0, compute the 2kp part of the self-energy in Matsubara frequencies and compare the prefactor for 
w n In |u;„| term with 1/2 of that in Eq. (|3.18|l . 
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FIG. 9: Hartree contribution to the self-energy for scattering at a Friedel oscillation. 



For a contact interaction, the second-order self-energy is 

H{k 1 u n ) = -U 2 J J ^ m G (k + q, u n + Sim) H P h(q, On). (CI) 

Assuming q — 2k F + q, where q is small, we expand tk+q as £k+q = — £fc + v F q + 2v F kp(l + cos6>), where # is the 
angle between k and q. As we already discussed in Appendix iBl only 9 near 9 = tt matter (i.e., typical q is nearly 
antiparallel to k), hence we can further approximate ek +q as 

ek+q ~ -efe + v F q + v F k F 9 2 , (C2) 

where 9 = tt — 9. Substituting l|C2|) into l|Blj) . we obtain, setting = 0, 

277 2 fc poo poo ^ 

SK) = 7 — f / dqdsi m / — — ^_^n ph fen ro ), (C3) 

where n p ^(q, r2 m ) is given by l|2.4(l . 

As an exercise, consider first a model case where n p ^(q, f2 m ) is static. To ensure convergence, we assume that 
the static behavior holds for Q, m <C f^O) where JIq is some ultraviolet cutoff (of order bandwidth), and for larger 
^m, n p /j(g, O m ) rapidly falls off. The angular integration in i|Cl(l reduces the range of integration over Sl m to 
— oj n < f2 TO < w n , hence at the smallest w„, E a w„. This accounts for the conventional mass renormalization. We 
now show that there are no non-analytic corrections to £ in this model. A static H p h(q, 0) is non-analytic only for 
q > 0, where H p h{q, f2 m ) = (m/2w)(l — (q/kp) 1 ' 2 ). Substituting non-analytic part of the polarization bubble into 1|B5|) 



introducing 9 = J r/v F k F cos <j), yfq — y/r/vp sin<^, and integrating over </>, we obtain for a potentially non-analytic 
part of the self-energy 

em r r <C4) 

One can easily make sure that this integral yields a regular ui term (determined by high-energy states), but no universal 
<x> 2 lnw-term. This implies, as we mentioned several times in the text, that static Tl p h(q, 0) does not give rise to a 
non-analyticity in the fcrmionic self-energy. 

It is instructive to distinguish this case from the impurity problem. If one of the interaction lines in Fig. la is replaced 
by an impurity line, as shown in Fig. [U] the diagram in Fig. la transforms into the Hartree diagram describing the 
scattering of fermions by Friedel oscillations produced by impurities. In the ballistic limit, \u> n \ r 3> 1, it suffices to 
keep only a single impurity line connecting G and Tl p h and also neglect disorder in G. For delta-correlated disorder 
with amplitude V, the analytic expression for the diagram in Fig. Stakes the form 

E(k,u n ) = -2UV Jj^ G (k + q,w„) Hp h (q,0). (C5) 

The particle-hole bubble is still static, but in distinction to IC1I we no longer have to perform a summation over 
frequencies. The non-analytic piece in E(w) is then given by, instead of Ea. (|B6|) . 

mUV f°° rdr 



=W— en; / —■ (06) 
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Due to the absence of the integral over fi OT , (|C4|) does yield a universal contribution T,(u> n ) oc — iui n \n(— iui n ) which 
comes from the lower limit of the integral over r. Upon analytic continuation, one obtains T,' R oc win \uj\ and T, R (lu) oc 
\co\ . The linear in uj form of T, R (co) is related to the Hartree part of the linear-in-T term in the conductivity at finite 
T [13 . 

We now come back to the electron-electron interaction, when a non-analytic- in-w„ behavior of E (w„) can be 
obtained if the f2 m — dependence is retained in H p h(q, Q m )- As with any logarithmic singularity, typical q should well 
exceed uj n /v F . We will see that typical O m are of order u n . Typical values of v F q well exceed then typical values of 
Q, m , and one can expand H p h{q 7 O m ) in powers of fl m /v F q. For q > 0, the frequency expansion of Tl p h(q, £l m ) starts 
at a constant and holds in even powers of fl m /v F q. We have already verified that the constant term does not give 
rise to an w 2 lnw-piece in E(w). At q < 0, however, the leading expansion term has the same \ fl m \ non-analyticity 
as the polarization operator near q = 0. The non-analytic behavior in frequency is crucial as it prevents one from 
eliminating a low-energy non-analyticity by closing the integration contour in the integral over fl m over a distant 
semi-circle in a half-plane where the denominator in (|B5J) has no poles. 

Expanding Tl p h{q, Q m ) at q < and fl m <C v F \q\ , we find 

u ph (q, n m ) = £ ( i - '"7' ) • (C7) 

Substituting this result into Ea. ljB5|l and keeping only potentially non-analytic piece, we obtain 

(2tt)4 J_ x m J_ x J v F q + v F k F e 2 -i(oJ n + n m ) v F (k F \q\y/ 2 
Introducing x 2 = —v F q and y 2 = v F k F 6 2 , we obtain from (|C8|I 

mU2 f°° jo m , f°° f°° dx d y 



EM - "3^ L«- |nJ /. I f-^-iktm - (C9) 

Introducing further y = y / tcos^/2, x = y / rsin</i/2 and integrating over <f> first, we obtain 

mU 2 f°° ,„ ,„ , f w dr 



8ir 4 v F J_ 00 J r J cos0- i(u + fi)/r 

= _i 8^| Sgn(w " + " m) A (r' + K + W/' " (C10) 

Evaluating the integral over r with logarithmic accuracy and integrating finally over f2 m , we obtain 

^H = - I ^->^ r - (en) 

This coincides with the half of Ea. H3.18(l for = 0. 

To further clarify this issue, we redo the calculation in a different way. Namely, we use the fact that for Q = — 2k+Q', 
and Q' small, the non- analytic part of the bubble H p h(Q' , fim) comes from the region of small Q" in the following 
integral: 

U ph (Q',Q. m ) = -J J 7^3F^ Gk+Q'V„ G -k+Q'+Q", w „+n m - ( C12 ) 

Now, we want to re-express the 2k F contribution as an effective Q = contribution. To do this, we substitute l|C12|) 
into (|B1|) and change the order of the integrations over Q' and Q . The non-analytic "2fcj?" piece in the self-energy 
then becomes 



) f" 2 / / d Y<!?y m G k +Q ',.„ + o m n(0'^ m ), (Ci3) 



where the effective particle hole-bubble 



tl(Q",Q m ) = ~ I I ® 2n y n G-k+Q",a> n G_ k+QI+Q » tUn+Um . (C14) 
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This II is a part of the particle-hole polarization bubble at small momentum transfer, which comes from the integration 
over small Q . We now show that for Q m <C vpQ , i.e., in the momentum/frequency range which yields the logarithm 
in the self-energy, the non-analytic part of IT(Q, fi m ) is a half of that in Ii(Q 1 O m ). This would again imply that the 
2fcp contribution to the self-energy coincides with the (non-divergent) E2 part of u q — 0" contribution. 

The calculation proceeds as follows. We set = and write £-k+Q' = — a;cos#i + 7a; 2 where x = vpQ', 7 = 
(2mvp)^ 1 , and 8\ is the angle between k and Q'. Similarly, e_k+Q' = — xcos6*i — y cos82+j(x 2 +y 2 + 2xycos(8i — 82)), 
where y — vpQ, and 62 is the angle between k and Q. As we said, we need to evaluate IT for 82 close to ±7r/2, and 
small y. We therefore neglect y 2 terms and set 82 ~ vr/2 for dcfiniteness. We assume and then verify that U/vfQ 
term in the polarization operator comes from 8\ near ±717 2 and linearize cos#i near these points. The integration 
over 8 1 is then straightforward, and performing it we obtain that the integration over u> n is confined to — Q m < u) n < 
(for definiteness we assumed that Q m > 0). The result is 

in m r . ( 1 1 



IL(Q,Q m ) = - / dp — + — — — , (C15) 

4ir 2 v F jy J \ cos 8 2 - 2p - iil m cos 8 2 + 2p - iil m J 

where we introduced p = jx. The integration over p is straightforward, and for small fi m and cos 82 the integral over 
dp yields in/2. Substituting this into (|C15(I we obtain 

n(Q,n m ) = if % (cie) 

It is essential that the momentum integral is confined to small p = Q' /kp (typical p ~ cos^2 ~ ^m/vpQ), and hence 
we are really restricting our momentum integral to small Q' . Comparing (|C16(1 and l|2.2|l we see that, as we expected, 
(|C16(l is a half of a non-analytic part of H(Q, f2 m ) at fl m -C vfQ. Another half obviously comes from the region of 
large Q' , which cannot be re-expressed as a ll 2kp contribution." 

2. An alternative computation of the self-energy, via U pp (q, fi) 

We discussed in the text that the second-order self-energy can be equivalently presented as a convolution of the 
fermionic Green's function and the particle-particle bubble 

d 2 qdil r 



S(w„) - -U 2 I I ; ™ Goi-kph + q, -uj n + Q m ) Tl pp (q, Q m ), (C17) 



where U pp (q 7 fl m ) = (m/2n) In B/ [\Cl m \ + y^m + (vFq) 2 
ing e-k F ii+q as —vpqcosd, we obtain for £fe = 



Substituting this H pp into the self-energy and expand- 



E(« n ) = r~o" / dQ m / d8 / dx In , (C18) 

8tt 4 w|, j_ 00 j j x cos e + i(n m -ui n ) \n m \ + v / n m + x 2 

Assuming, as before, that typical fi TO are of order u) n , while typical x = vpq are much larger, we can further expand 
under the logarithm and obtain 

tj2 rOC />7T />M / -I 

mil /,„,„,/,„/, 1 



S K) = / dn m \n m \ / dd / dx — — r. (ci9) 

8tt 4 w|, J_ oq 7o 7o x cos 6 + i(n m -ui n ) 



The integration over 8 yields 



mil 



2 r°° r w 



E(u) n ) =i »' / dVL m \Vt m \ sgn(fi m - uj n ) I dx (C20) 

Evaluating the integral over a; to logarithmic accuracy, we finally obtain 

mil 2 W 2 

S(u ra ) = -i — 2 u) n \n . (C21) 

This coincides precisely with Eq. I|B11|I . 
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APPENDIX D: EVALUATION OF T.' R (uj,e k ) ON THE MASS SHELL 

In this Appendix, we present the calculation of the real part of the fermionic self-energy on the mass shell. We 
will be only interested in the non-analytic piece of the self-energy. The non-analytic part of Y,' r (uj) is simply twice of 
Y,' 2 (tu), where, according to Eq. (|3.37|) . S^w) can be written as 



(Dl) 



where 



z(lu,t) = [ dnn v [ 

J -oo J 



dE 



UJ 



E 2 -uj 2 
2VL + E-uj 



2VI + E + UJ 



hi 



, n , n+E 

coth — — — tanh — — — 
2T 2T 



\{2VL + E) 2 



W 2 



(D2) 



We first find Z(iS) at T = 0. The term with coth and tanh functions restricts the integration over 17 to the interval 
— E < fi < 0. Introducing the rescaled variables E — ujz and Q = —luzx and assuming for definiteness that u> > 
(and thus z > 0), we obtain 



Z(lu) = 2u> 



dz 



xdx 



z In 



z{2x-l) + l 



z{2x - 1) - 1 



+ \n[z 2 (2x- l) 2 - 1] 



(D3) 



Introducing a new variable via y = 2x — 1 and eliminating terms that vanish by parity we obtain, instead of Eq. (|D3|> , 



Z(uj) = LO 



dz 



dy 



zy In 



zy+1 



zy-1 



ln(z 2 y 2 



(D4) 



The integration over y is now straightforward, and performing it we obtain 



Z(uj) = uj 



dz 



1 \z 



In 



z+l 



z- 1 



+ lnz 2 - 1 



dz 

Yz 



In 



z + l 



z - 1 



(D5) 



Finally, we use the values of the following integrals 



f 

Jo 



dz 



lnz 2 



1 



7T~ 



dz 1 



1 Z 



In 



z + l 




IT 2 


Jo 2z 


z + 


1 


z-1 




~ 4 ' 


z — 


1 



7T~ 

T 



(D6) 



Substituting these results into Eq. l|D5(l we obtain 



Z(lo) = uj- 



(D7) 



Substituting this further into Ea. lDl|) we reproduce Eq. 

We next consider finite T. As a first step, we show that one can safely replace cothi7/(2T) by tanhfi/ (2T) in (|D2 
Indeed, this replacement changes Z(uj) by 



n 

rffi — n r 

, sinh M 



dE 



E 



hi 



2fi + E - uj 



20 + E + uj 



In 



|(20 + £Q 2 -uj 2 \ 
W 2 



(D8) 



The integration over E in Z ex t ra (w,T) is straightforward and performing it we obtain 



Z e xtra,(^J, T) 



dn 



il \n 2 \2fl + uj\ - In 2 \2Q - uj\ 
sinh p w 



(D9) 



This integral obviously vanishes as the integrand is odd in fi. 
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Next, one can readily check that in the expression for Z, obtained by replacing cothil/(2T) — > tanhO/(2T), i.e., in 

dE 



Z(lo,T) 



dun v 



10 



E 2 

20, + E - uj 



2H + E + UJ 



In 



(tanh^-tanh^ 
V 2T 2T 



|(2fi + Ef -, 
W 72 



(D10) 



the integrand vanishes at large |f2|, E. Hence the integration can be performed in the infinite limits and Ea. (|D10|) 
can be rewritten as a difference of two terms with the same argument of tanh, upon changing in the second term to a 
new variable f2 + E. Carrying out this procedure, introducing new variables, and converting the f2 integration to the 
integral over positive f2, we obtain 



2fl \ 



where 



*(a)=V 



dxx 



1 



f°° fl 
Z(to,T) = / dfttanh — * 
Jo 2T 

a 2 - (x - I) 2 , (a - l) 2 - x 
~) In 



• In 



{x + lf 



{a+lf-x 2 



In 



x) 2 - 1 



{a + xf 



(Dll) 



(D12) 



The integration over x is tedious but straightforward, and yields 

*(x) 



7T 2 §, for a < 2; 
7r 2 , for a > 2. 

Substituting this into Ea. (|Dll|) and integrating over Vt we obtain 

7T 2 |w| 



Z(w) = A + 



-9 



T 



(D13) 



where A < is a (formally infinite) constant which is irrelevant to us as it accounts for the high energy contribution 
to a linear in u> term in T,' 2 (lo), g(x) is universal scaling function 



1 



12 



Li 2 



(D14) 



and Li2(x) is a polylogarithmic function. This is the result we cited in Eg . Q3 .42(1 . 

At x = oo, i.e., at T = 0, we have 5(00) = 1 and thus Z{u>) = (it 2 /2)\to\. This coincides with Eq.|D7|). In the 
opposite limit of \u>\ <C T, we use property 



Li 2 (-e-) = E 
fe=i 



(~e- a: ) A 
/fc 2 



12 



xln2 + C(a; 2 



Substituting this into Eas. (|D14fl and (|D13|I wc obtain that up to a constant, 

Z(u < T) w 2tt 2 In 2 T 
Substituting this further into Ea. ljDll) we obtain 



m[/ 2 In 2 
8tt 2 4 



(D15) 



(D16) 



(D17) 



This is the result we cited in Ea. H3.39|l . 

As an independent verification, we reproduced i|D17|) by computing the temperature derivative of Z(u>) in the limit 
lu —>■ 0. [It is essential to take the limit, not just set to = 0]. Evaluating the derivative, setting to — > 0, introducing 
dimensionless variables, and eliminating the terms which vanish by parity, we obtain 



dZ(w,T) 
df 



= 4 



dxx 
cosh 2 : 



dy 

y 



in 



y 



1 



y-i 



The integral over x gives In 2, whereas that over y yields, upon integrating by parts, 

dy In y ~ 



V 



1 



7T~ 

T 



(Dig) 



(D19) 



Combining the two terms we obtain dz< p^ T ^ 
with Ea.0T5l. 



2-7T 2 In 2, i.e., up to a constant Z(u <C T) = 27r 2 In 2 T. This coincides 
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APPENDIX E: 2k F CONTRIBUTIONS TO DIAGRAMS 1 AND 3 IN FIG. 5 

In this Appendix we present explicit calculations of the 2k F — contributions to diagrams 1 and 3 in Fig. [3] 

1. 2kF part of diagram 1 

We first verify that the non-analytic 0(|Q|) term that results from the 2k F non-analyticity in the particle- hole 
bubble is indeed the same as the contribution from the q = non-analyticity. For 5x\~ (Q, 0) we obtained in Ij4.18jl 

A 9=o (n ^ 2 (mU\ 2 \Q\ 

Now we explicitly evaluate Sx\ kF (Q, 0). The general expression for the diagram 1 is 

/d 2 k d 2 a dtndQ 
7 Gg(k, c)G (k + Q, c)G (k + q, cj + n)U(q, fi). (E2) 

For q w 2fcp the quasiparticle energies can be approximated by 

€k=v F {k-k F ), ek+Q = £fc + vpQcosOi, ek+ q = -£k + v F q + 2v F k F (l + cos0 2 ), (E3) 

where q = q — 2k F , and 9\ and #2 are the angles between k and Q and between k and q, respectively. As we have 
said several times before, the 2k F non-analyticity comes from internal fermionic momenta in the particle- hole bubble 
that nearly coincide with the external one. In our notations, this implies that #2 is close to tt. We can then expand 
in cos 6*2 upon which ek+q reduces to tk+q = — £fe + v F q + v F k F (jr — Q2) 2 ■ Substituting this expansion into (|4.9|1 . 
integrating over tk and then over to (this requires more care than for the q = case), and introducing dimensionless 
variables q = q/\Q\, ui = Cl/(v F \Q\), k F (ir — Q2) 2 = \Q\Q 2 and polar coordinates as q = rcos<j>,uj — rsin</>, we obtain 
from (|4.9|l 

jxf , {m _ mm» r mw * r rdr r * (E4) 

7T 4 W F J J J Q COS^fc*! 







cos6»i , 6 2 +re lc t> + cos6»i 
In 



The polarization operator is now given by l|2.4|) which in the new variables takes the form 



(E5) 




11(0)=^-- 1- cos^ . (E6) 



Performing the integration over r and keeping only the contribution which comes from low energies, we again find 
that only the non- analytic piece in 11(0) contributes to order \Q\, and this universal contribution is 

3ir b v F J 2 J (6» 2 + e l <t>Y 

r 9 2 + 

x / d0! cos6»i In . (E8) 

Jo cos 01 

The integral over 9\ yields wr. Evaluating then the integral over 9 2 , we obtain 

r 2fc P ^ ^ m 2 U 2 \Q\ r 4> . 50 

5?& kF (Q,0) = - 4 m / d0cos^ sm^ 

2 (mU\ 2 \Q\ _m 

Comparing this result with Eq. (|4.18|) . we see that the two expressions are indeed equal. We emphasize again that 
in order to obtain this result, one has to include the frequency dependence of IL(q,u>) near q = 2k F . Had we replaced 
II(<7, cu) by its static value (II(q, 0)), we would have not obtained Eq. (|E9|) . 
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2. 2/cf part of the diagram 3 



In an explicit form, 



SxsiQ, 0) = -4[/ 2 J l|-^^ G (k, w )G (k + Q, uj)G (k + q, u + fl)G (k + q + Q, u + fl)U ph (q, fl), (E10) 

Assuming that q is close to 2k F and expanding quasiparticle energies as in Ea. (|E3|) we obtain after rescaling the 
variables and restricting with only the non-analytic part 



S X3 (Q,0) = xo ^-—J_j x j o dfl 

1 



( \fx + ifl + \/x — iflj 



oo poo 



d6 dy dz du>-. 
Jo J-oo J-oo (z - iu)(z + cos 6 -ilu) 

1 

(z — x — y 2 + i(oj + fl)) (z — x — y 2 + cos + i(oj + fl)) ' 
where xo = tu/tt. Performing the integration over z first we obtain after straightforward manipulations 

m 2 U 2 \Q\ 



(Ell) 
(E12) 



SxT F (Q,o) = -xo : 



d6 I dylra 
' Jo 



, dx I dfl 



(^Jx + ifl + Vx~-iflj 



o ' (x + y 2 - i{uj + fl)) ((x + y 2 - i(w + fl)) 2 - cos 2 



(E13) 



Introducing then q — rcostp and fl = rsin7r such that (V x + ifl + yx — ifl) = 2^/r cos cf>2 and rescaling ui — ► rw, 
y — > \/ry, we obtain 



7T 5 /CF Jo 



Illl 



oo poo 



1 



m ^/o Jo r Jo d9 (e-^ + y 2 - iLu) {r 2 {e- l <t> + y 2 -iu) 2 -cos 2 < 



(E14) 



Introducing then p — r(e + y 2 — iui), replacing the integration over r by the integration over p, and restricting with 
the universal contribution from the lower limit of the p— integral, we obtain, after integrating over p and then over 9 



Sx? F (Q,0) = -2 Xo - 



2 u 2 \q\ r 



7T 4 k 



"J f--K /> OO pOO 

— / d4>cos4>/2 / dy dujRe 

F Jo Jo Jo 



1 



(uo + i(y 2 + e-^)) 4 



(E15) 



The integration over u> is now straightforward. Performing it and then evaluating the integral over y we finally obtain 



Sxi kF (Q,o) = xo 



m 2 U 2 \Q\ r 



in 3 k F Jo 



<j> 2 fmU 

, COS - S1 n-= Xo - — 



(E16) 



This is the result that we cited in the text. 



APPENDIX F: 2k F CONTRIBUTION TO Xs (Q = 0,T) FOR A STATIC LINDHARD FUNCTION 



In this appendix we show that the thermal smearing of the static Lindhard function by itself does give rise to 
a linear-in-T term in the uniform spin susceptibility, but does not account for the full linear-in-T dependence of 
X s (0,T)-the latter also contains a contribution from finite frequencies. 

The computation proceeds as follows. Because a static polarization operator can be viewed as an effective interac- 
tion, diagram 1 can be re-expressed as the first-order self-energy insertion (see Fig lTU)) 



§Xi,static = -4T^ J 



d k 3 

^[G(k,U) n )] S e ff(efc), 



(2tt 



(Fl) 
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FIG. 10: Diagram 1 of Fig. 3 as the first-order self-energy insertion. 



where the effective self-energy is given by 

E cff (e fe ) = 2[/ 2 rV f ^n(?,0,r) G (k + q, W „) = 2U 2 [ ^n(,,0,T)« f (ek +q ) 
„ J (27r) J (2tt) 

This self-energy is obviously independent of u>„,. Although the static polarization operator H(q,0, T) is not known 
exactly, it can be cast into an intergal form [31| convenient for further calculations. We have 



k 2 F f^ 2 - 1 dz f, 1 + z \ 1/2 ' 



SmTj^ C osh 2 ^ V (q/U F ) 



AiuT 



1 - 



(F2) 



Re-writing [G] 3 = (1/2) d 2 G/de 2 , summing over u n with the help of an identity TJ^^ G(k,u> n ) = n F (e^) — 1/2 
where n F (z) = (e z / T + l) is the Fermi distribution function, and integrating by parts twice, we obtain 

SXhatatic = -Xo D [ de k n F (e fe ) - Scff J £fc ^ (F3) 

J-oo ae k 

where Xo D — fn/it. A non-analytic temperature dependence of 8\i, static is due the region of q near 2k F , where 
Tl(q, 0, T) is singular. Expanding, as before, £k+q near q = 2k F and along the direction of q nearly antiparallel to k 
because only these q contribute to the non-analyticity, we obtain 

e k+q = -ejt + v F (q - 2k F ) + v F k F (it - 6) 2 , (F4) 



where is the angle between q and k. Substituting £k+q into l|F2(l and rescaling variables, we obtain for the effective 
self-energy 

Tj2 u2 / 07^ \ 2 l-oo j-oo r x I \l/2 

S off (6 fe ) = -^-^ / dx dy dz [X g n F (-e k + 4:T(x + y 2 )). (F5) 

27T 3 \&F J J-oo Jo J-oo COSh Z 

Substituting this self-energy into l|F3[) . evaluating the derivative and further rescaling variables we obtain 

5 XhstaUc = ~xl D f^V ^rZ, (F6) 



4tt J E 



F 



where 



and 



Z=-\ da I dbT^ ^(fe-c) 

- f ' , io ^ cosh 3 (6-c) V ; ' V ; 



•'<»■'"=/ e 4(nJ-,) + 1 - ( F8 ) 
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The last integral can be easily evaluated and yields 

= J£rrr (F9) 

Substituting this result into (|F7|I . introducing c = y/c, and b = a 2 + b and integrating over c and a using polar 
coordinates, we obtain after straightforward calculations 

j.-.r *- V (Flo 

./_^ e 46 - 1 cosh 2 6 



Carrying out the last integration, we finally obtain Z = — (1 + 7r 2 /4) and 

^*uu c = xl D f^V - fl + V) 'FID 



47r / ej? 

Comparing this result with our 5x\~ 2kF = ^Xi _ ° = (l/2)^Xs(0, J 1 ), given by Eq. H4.54fl . we see that they differ in 
that Z 7^ 1 . This discrepancy shows that the frequency dependence of the polarization bubble does contribute to the 
non-analytic piece in the thermal static uniform susceptibility. 
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